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1. Comparison of Plasma and Liquid Turbulence 


AS AN introduction to the whole problem of plasma turbulence, we start 
with a comparison of plasma turbulence and the turbulence in liquids. The 
turbulence of liquids has been the subject of an enormous number of 
investigations during the last 30 years, and therefore the basic concepts 
are well known. The field of plasma turbulence has been developed 
only in the last decade both experimentally and theoretically. Never- 
theless, it is surprising to see how much research is now devoted to plasma- 
turbulence problems as compared to liquid-turbulence problems. On the 
one hand, this is due to a very large activity in the experimental and theo- 
retical study of plasmas in recent years; these activities are partly due 
to the applications and explanations of the phenomena found in labora- 
tory investigations of plasmas in high electric and magnetic fields and 
of shock waves, as well as those phenomena found in plasma heating 
or in geophysical and astrophysical applications. On the other hand, 
there exists a physical reason which allows us to give a complete theo- 
retical description of the most important turbulent motions connected 
with so-called plasma oscillations. This is the existence of some kind of 
elasticity in collective motions of plasmas, which is absent in incompres- 
sible liquids. For example, if in a plasma, a sheet of plasma electrons is 
displaced over a distance d (see Fig. 1) the charge separation provides 


t+tt¢te+44 





Fic. 1. 
Plasma oscillation arising from a charge sheet separation. 
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a force which tends to prevent the charge separation and oscillations 
around the neutrality position arise. These are the so-called Langmuir 
plasma oscillations which have a frequency 


ay wy 


where n is the mean density of the electrons and m, their mass. In the 
presence of an external magnetic field a frequency of the order of the 
gyrofrequency of electrons 


One = -->5 (1.2) 


also occurs; other frequencies, involving the ion parameters, are also 
possible (Stix, 1962; Ginzburg, 1970). 

In incompressible liquids the eddies have no special frequency and 
their frequency is determined by their interaction with the other eddies. 
Turbulence is usually connected with the non-linear interaction of 
collective motions. Owing to the elasticity of plasma motions there 
occurs a small parameter, which is the ratio of a period of oscillation 
1/@ to the characteristic time t of the non-linear interactions, that is, 
we have 

oo 1. (1.3) 

Condition (1.3) is called the condition for weak turbulence. The theory 
of plasma turbulence was completely developed for the case when (1.3) 
holds, and this condition seems to be valid in most experimental investiga- 
tions of plasma turbulence (Kadomtsev, 1964). 

In Fig. 2 we show the possible branches of collective plasma motions. 
Plasmas can also have collective motions similar to those which incom- 
pressible liquids have, namely eddies. These motions exist for @ « » (v is 
the two-particle collision frequency) and for w/k « v,, where v, is the 
sound velocity. This region is the very small shaded region in Fig. 2. 
Most of the collective motions shown in Fig. 2 are collisionless. Let us 
mention that c,./»- is usually a very big number (as a definition of the 
plasma as a state of matter implies) from 104 in laser produced plasmas 
to 108 in most laboratory gas discharge plasmas, and to 101-10! in 
astrophysical conditions. We see, therefore, that the concepts of plasma 
turbulence contain a generalisation of the concepts that were used earlier 
in the case of liquid turbulence, 
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Fic. 2. 
Possible branches of collective plasma motions. M and MS indicate the 
fast (w = kv,) and slow (w = kv,) magnetohydrodynamic waves, A Alf- 
vén waves (w = k,v,), W whistlers (w = wg,k2c?/w?,), S ion-sound waves 
(w = kv,), i plasma ion oscillations, h hybrid plasma oscillations (@ = 
w,-k,/k), | Langmuir plasma oscillations (w © w,.), and t, andt, ordi- 
nary and extra-ordinary transverse waves. Here v,>>v,, where v, is the 
Alfvén velocity (v, = H/+/(41nm;)), v, the sound velocity (v.= /(T./m,)); 
pe >> Waa; and the external magnetic field is assumed to be in the z- 
direction. 


The other essential difference between plasma and liquid turbulence 
is the presence of electric and magnetic fields in collective motions. 
Therefore, if in a liquid the stochastic variables are the mean particle 
velocity of particles » or the density gin plasmas the electric fields also 
become stochastic properties. 

In the general case one can divide the observables into two parts—a 
regular and a stochastic part 


vo = ofe8+- pstoch 5 (1.4) 
E —_ Eres 4 Estoch . (1.5) 
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By definition 
CE slosh) = 0, (pstoch), = 0, (1.6) 
where brackets indicate averages (for example, over a statistical ensemble). 
As in liquids the statistical properties of the stochastic fields are given 
by correlation functions that are measured in most experiments of 
plasma turbulence: 


i ((Esoch(r, t)-Esoh(r’, 1) = | dea dK Ty, e+ HED, (1.7) 


Equation (1.7) is written for the case of a stationary and homogeneous 
turbulence, and when r = r’ and t = f¢’ the quantity J,,,, gives the fre- 
quency and k-dependence of the energy of the electric field of turbulent 
motions. 

Usually as, for example, for all turbulent motions shown in Fig. 2, 
the whole energy W,,., of turbulent motion, which consists of the energy 
of the turbulent electric fields and the energy of the particle motion, is 
proportional to J, .,: 

Wr,w => Ok, wlk, w+ (1.8) 

Equation (1.8) is an approximate one. It is approximately valid, if (1.3) 
is valid. The w-width of W,,,,, that is, 4w, describes the characteristic 
time of the correlations: 


1 
Tcorr ~ Ae (1.9) 


The correlation curve W,,,, of a collective plasma motion as a function 
of for fixed k has a maximum near the frequency of elastic response 


(Fig. 3). 





Oy 
Fic. 3. 
Sketch of the correlation curve in a turbulent plasma. 
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The half-width Aw of this curve is much less than w(k), if the condition 
(1.3) for weak turbulence is fulfilled. As in liquids the integral of W,,., 
over w is called the turbulent spectrum 


W,.= i Wr, wo do. (1.10) 


For the case of isotropic turbulence it is useful to introduce a quantity 
W,, normalised to the modulus of & (|K| = k): 


Wi = 4nk?W,,, C1 


[weak = w. (1.12) 
0 


For magnetic types of oscillations it is also useful to introduce a quantity 
W.., 2, which satisfies the equation 


[Pao [ a22.0 ex (1.13) 
0 


where 9 is the solid angle. 

Now let us come to the problem of the spectrum of turbulent liquids. 
Suppose, as a very simple example, that we have a liquid flow in a pipe 
of dimension Lo (Fig. 4). 





Fic. 4. 


Sketch of the velocity distribution in a liquid flowing through a pipe. The 
solid line corresponds to laminar flow and the dashed line to turbulent 
flow. 


If the mean velocity (or the Reynolds number) of the liquid is low 
enough, there exists a laminar flow which forms a velocity profile such 
as shown on Fig. 4. If the velocity reaches some critical value, such 
a profile becomes unstable and eddies of dimension Lo are excited. These 
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eddies become also unstable, excite new ones with smaller dimensions, 
and so on. In the developed turbulent regime, eddies with dimensions 
less than Lo are present, so that one can ask for the distribution of the 
turbulent energy over the dimensions of the eddies, that is, the energy 
that is carried by eddies with dimensions in the interval between / and 
1+<dl. Instead of / one can introduce the wave-number k, 


k=" (1.14) 


and ask for the distribution over k. This is the same as that given by 
equation (1.12). The region k ~ ko = 22/Lo is usually called the energy- 
containing region. The diminishing of the eddy dimensions described 
above is usually considered as a result of their non-linear interactions 
and leads to an energy flow in k-space from k = ko to higher k. For 
very large k viscosity becomes important and the turbulent energy of 
the eddies is decreased by viscosity (see, for instance, Landau and 
Lifshitz, 1959). 

We show in Fig. 5 the turbulence spectrum of incompressible liquids. 
This spectrum is stationary. That means that for each k there is an energy 


Wy 





Fic. 5. 
The turbulence spectrum of incompressible liquids; Q is the energy input, 
W,, ~ k~5!3 the universal Kolmogorov law. 
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balance between the eddies which enter and leave the interval dk. In the 
energy containing region (k of the order of ko) the balance is maintained 
by the excitation of eddies and their transformation to higher k by non- 
linear interactions. In the middle part of the spectrum, called the inertia 
region, the balance is only caused by the non-linear interactions. That 
means that the energy which enters the interval dk from lower k-values 
is compensated by energy outflow from dk to the higher k. In this region 
the spectrum is universal; it was first found by Kolmogorov (1941) to 
vary as ~ k—5/3 by means of dimensional analysis and is known as 
the Kolmogorov spectrum. In the region of high k, when the damping 
of the eddies becomes comparable with the energy input due to the 
energy inflow from small k-values, there is a balance between energy 
inflow and damping. If one raises the energy input Q the energy spectrum 
does not change (see Fig. 5, Q2 > Q). All the characteristic elements of the 
formation of the stationary turbulence shown above for liquids are also 
found in the case of plasma turbulence (Tsytovich, 1969a). These are 


1. The excitation of collective motions in some interval in k-space. 

2. The non-linear energy transfer from the excitation region to the 
region where the turbulent oscillations are damped. 

3. The damping of the oscillations. 


The difference lies, of course, in the actual mechanisms of excitation, 
energy transfer, and damping. 

Since damping is the final result of the history of any turbulent oscilla- 
tion excited ina plasma, it is useful to classify the possible types of plasma 
turbulence by mentioning the most important damping mechanisms. 
There are new mechanisms of damping that have no analogy with the 
damping of eddies in liquids, and, therefore, in a plasma new types of 
turbulent motions can exist. First, it should be mentioned that the 
turbulent oscillations can be damped by binary collisions, and this kind 
of damping is in some sense analogous to the damping in liquids. (Al- 
though in a collisionless plasma w >> v and in liquids w< v.) This kind 
of damping can be essential, for example, in a partly ionised plasma, or 
in the case when the turbulent oscillations are transformed by non- 
linearity from the regions of intensive collective damping. Such a type 
of turbulence—turbulence dissipated by binary collisions—can have a 
large energy stored in the turbulent motions because the damping is 
usually small. 

Secondly, in plasmas there exists a new type of damping—the collective 
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damping. An example of such damping is the well-known Landau 
damping. The presence of such damping is connected with the electric 
field that exists in turbulent plasma motions and, therefore, with the 
possibility of the interaction of such a field with the charged thermal 
particles of the plasma. It is essential also that one can associate with the 
turbulent motions definite phase velocities 


ia (1.15) 


Because of correlation broadening this correspondence is only approxi- 
mate (see Fig. 3) but the uncertainly is of the order of the small parameter 
(1.3). If then v7, becomes of the order of the mean thermal particle 
velocity, all thermal particles are resonant with the waves, that is, it can 
gain energy from the wave and, therefore, the turbulent oscillations are 
heavily damped (Landau, 1946). 

In a magnetic field there can also exist cyclotron resonance, when the 
frequency of the field of turbulent oscillations in the reference frame 
where the particles have zero velocity along the magnetic field is equal 
to the gyration frequency of the particles or to its harmonics, This is the 
so-called cyclotron resonance, or cyclotron damping. This kind of 
damping can exist near the electron-cyclotron or near the ion-cyclotron 
frequencies. 

If the energy of turbulent motions flows in the direction where a region 
of collective dissipation exists, the heating of plasma is very essential. 
Such a kind of heating is usually called turbulent heating. From the 
considerations given above, it is obvious that the rate of such heating 
depends on the rate of energy transfer and therefore depends on W, or 
the energy input. 

Thirdly, in a plasma there can exist also a new type of damping due 
to acceleration of charged particles. Such an acceleration is a stochastic 
one and is also deeply connected with the electromagnetic nature of the 
plasma oscillations. One can call this also a heating of a small fraction 
of the particles which are resonant with the oscillations, that is, with 
the maximum in Fig. 3. The thermal particle can be non-resonant, and 
then the oscillations are not heavily damped. Because of the long life of 
the oscillations, fast particles can receive enough energy and as a result 
the oscillations can be damped. Such particles occur naturally in the tail 
of a Maxwellian distribution in a plasma. Their injection into an acceler- 
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ation regime (that is, one in which they are able to resonate with 
oscillations) can be due to the non-linear interactions of the oscillations 
or to a magnetic type of turbulent waves. 

Fourthly, there can also be a transformation of the turbulent energy 
into electromagnetic radiation. It is due to the electromagnetic nature of 
turbulent plasma motions and their non-linear interactions. This kind 
of turbulence is radiatively dissipative turbulence. 

The excitation of turbulence in plasmas can be due to a change of 
sign of the damping coefficient as a result of the anisotropy of particle 
distributions. Thus the same collective effects that give damping in the 
isotropic case can produce instabilities, if the particles are distributed 
anisotropically. The turbulence can be excited also by fast particles and 
the best known example of this is the plasma beam instability. Also it is 
possible to excite turbulent motion by radiation or by a high-frequency 
electromagnetic field. This has an analogy with the non-linear optics 
problem of excitation of supersonic waves in solids (Tsytovich, 1967). 

The existence of charged particles in a plasma is very essential for the 
nature of the non-linear energy transfer. In liquids the eddies interact 
only with one another or, in other words, the whole energy of the tur- 
bulent motion is conserved and is only transformed from the biggest eddy 
to the lowest, or from smaller to larger k. The entropy in this kind of 
process increases because the phase volume that is proportional to the 
volume in k-space increases. 

In the presence of charged particles and due to the electromagnetic 
nature of plasma oscillations there exists the possibility of interactions 
and energy exchange between the turbulent motion and the particles. 
One of the most important of such processes is the induced scattering 
of oscillations by particles. One such process is shown in Fig. 6. In such 
a process the particles can gain energy or, in other words, are heated. 


Fic. 6. 
Scattering of oscillations by particles. 


TSY 2 


10 THE THEORY OF PLASMA TURBULENCE 


The particle entropy increases because of this heating, and to compensate 
this increase the phase volume occupied by the waves can decrease. In 
other words, the energy of turbulent motions can be transformed from 
higher to lower k. This is the case, for example, for Langmuir plasma 
turbulence, the spectrum of which is shown in Fig. 7. 





Fic. 7. 
The spectrum of Langmuir turbulence (7, = 7)). 


The turbulent oscillations created at high k (where the source of the 
turbulent energy is located) are transformed to lower k. The spectrum 
~ k-52 is due to the energy balance when this transformation is pro- 
duced by induced scattering by electrons. This spectrum exists only for 
k > (m,./m)'®kp in a narrow region. The spectrum W, = const = k7’, 
v = 0, due to the induced scattering by ions is valid up to 


1/2 
oF (Fn aes. fre (1.16) 


The turbulent energy is transferred to smaller k up to k, (each step of 
such a transfer is of order k,, so that one can only roughly say that the 
spectrum is flat on the average in a k-interval larger than k,).In the region 
k<k, the whole energy has come to small k and the eddies visually 
speaking must have very rapid collisions with one another because the 
phase volume is small. This repulsion is described by the interaction of 
plasmons which is diagramatically shown in Fig. 8. 
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uy f, 
Fic. 8. 


Plasmon collisions; 1, and | describe the plasma state before and 1; and I’ 
the state after the collision. 


As the collisions of particles produce a Maxwellian distribution, the 
collisions of plasmons produce a maximum in Fig. 7. The essential 
difference is that there is continuous flow of plasmons from higher k into 
this region. The damping of plasmons in the region of the maximum in 
Fig. 7, which can be called an energy-containing region, can be due either 
to ordinary binary collisions or to radiative losses. In the second case, 
the Langmuir turbulence can be radiatively dissipative (if the number of 
fast particles accelerated is small). 

The spectrum in the asymptotic region, k,>> k>>ko, is ~ k~”, where 
2:8 < v < 4, depending on the energy input Q. If the input of the turbulent 
energy Q is raised the spectrum remains the same in all above-mentioned 
regions except the last where » tends to4,ifQ increases. The position of 
the maximum of the spectra, ko = 2%/Lo, comes down in k-space if the 
dissipation is due to ordinary collisions and is unchanged if the dissipation 
is due to electromagnetic radiation. This picture was found by several 
authors (Pikelner and Tsytovich, 1968; Liperovskii and Tsytovich, 1969) 
by an analytical solution of the complicated non-linear equations of 
weak turbulence and checked by numerical computations by Makhankov 
et al. (1970). 

We shall describe the properties of Langmuir turbulence in more 
detail later on. We mention here that the energy transfer in this case is 
the opposite of that in liquids, and that this can be possible only if the 
particles are heated. Such a heating is a stochastic process. One must 
distinguish this from the turbulent heating where the energy of turbulent 
oscillations is directly drawn into the dissipative region. For Langmuir 
turbulence we have no such heating in the absence of magnetic fields 
because the direction of energy transfer in the case of isotropically dis- 
tributed particles is such that the turbulent oscillations are drawn away 


2* 
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from the region of Landau damping. From the energy-conservation law 
it is possible to see that the energy of plasmons becomes thermal particle 
energy. Because the number of plasmons is conserved, the larger the 
difference between the maximum and the minimum frequencies of the 
branch to which the energy is transformed the larger the heating. There- 
fore, the rate of stochastic heating depends on this difference as well as 
on the time needed for the transfer of waves from maximum to minimum 
frequencies. As can be seen from Fig. 2, the Langmuir waves undergo 
a small change of frequency along the branch and the low-frequency 
waves such as whistler and ion-sound waves a large one. On the other 
hand, however, the interaction of whistlers is rather small compared 
with the interaction of ion-sound waves and thus the transfer of the 
energy across the same frequency difference needs a longer time. The 
Langmuir turbulence can transfer the energy more quickly than the ion- 
sound turbulence. Both Langmuir and ion-sound turbulence are appro- 
priate for stochastic heating. 

It is well known that for liquids in the turbulent regime the viscosity 
becomes anomalously large and the profile of a liquid flowing through 
a pipe is flattened. Similar effects are observed in turbulent plasmas. For 
example, there occurs an anomalous diffusion in a plasma confined by a 
magnetic field, an anomalous resistivity when an electric field is applied, 
and so on. One can, therefore, introduce the effective collisions due to 
such processes. The physical nature of such collisions in the weak- 
turbulence case is very obvious, the collision frequency is the inverse 
of the characteristic time of energy transfer along the turbulent spectra. 
Thus one can expect that this characteristic collision frequency is detect- 
able if a small-amplitude wave with a frequency smaller than that of the 
characteristic turbulent collision frequencies is pushed onto the plasma. 
In this case, the skin effect of the electromagnetic wave can be due to 
such turbulent collisions. Indeed, a growth of the skin-depth when 
turbulence is excited was found experimentally. 

Generally speaking, one can say that the electromagnetic properties 
of a turbulent plasma governing the penetration ofa small linear perturba- 
tion is quite different from those of a non-turbulent plasma if the fre- 
quency of the perturbation is much less than the turbulent collision 
frequency. A new type of wave can appear and also a new stability pattern, 
that is, the stabilisation of some waves which are unstable in the non- 
turbulent regime and the existence of new instabilities. As an example, 
we mention the possibility of the excitation of a magnetic type of per- 
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turbations (Tsytovich, 1968a) in pure Langmuir turbulence, that is, in 
the case when only electrostatic fields are present. This effect has an 
analogy with the well-known effect of the excitation of a magnetic field 
by the turbulent motion of conducting liquids—the Batchelor effect. The 
scale of the characteristic dimensions or wavelengths of the excited 
magnetic fields in the case of Langmuir turbulence is much smaller than 
in the case of the hydrodynamic dynamo. 


2. General Problems of the Theory of Plasma Turbulence 


IN THIs chapter we shall give a summary and general description of some 
problems of the theory of plasma turbulence. 


2.1. Excitation of the Turbulent State 


Before discussing turbulence itself, one should define the kinds of 
motion which can be considered turbulent and how the motions can be 
excited. Usually the turbulence is connected with a stochastic variation 
of some measurable quantity, for example the electric potential g. In 
experimental investigations the time-dependence of the potential fluctu- 
ations can be measured. One can find many examples of such measure- 
ments in which @ changes in a very complicated manner with time (see 
Fig. 9). 


<7 
t 
Fic. 9. 
A typical experimental variation of the electrical potential! in turbulent 
motion. 


Does this mean that the potential is a stochastic variable? Generally 
not! For @ to be a stochastic variable, the measurement should not be 
reproducible by repeating the same experiment under the same macro- 
scopic conditions (which must be carefully defined). When is stochastic 
in this sense one can ask: What is the cause of the irreproducibility ? 
The answer is that small changes in the initial conditions change the 
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behaviour of » considerably. This is well known in the statistical theory 
of molecular motion. Plasma-turbulence theory deals with the collective 
motions of a plasma, such as plasma oscillations. There are no difficulties 
in describing this kind of motion such as there are in describing the 
individual particle motions and general statistical concepts can be 
applied to describe the turbulence of plasmas. 

For plasma oscillations, which in zeroth approximation could be con- 
sidered as propagating waves, the stochastic behaviour lies in the ran- 
domness of their phases. Indeed, the initial value to be specified to match 
the potential of such a wave is its phase. Randomness of phase means 


(p) =0. (2.1) 


The average is over the statistical ensemble of the expcriments mentioned 
above. In the general case, there is an applied or regular component of 
gy as well: 

Q= gree + gpstoch 7 (pstoch, = 0. (2.2) 


For experimentalists it is necessary to mention that the theory predicts 
the statistical average value which is to be found by averaging the results 
from many equivalent experiments. 

The question now is, why do small changes in the initial conditions 
change so appreciably the time behaviour of the measured quantity? The 
answer is that such initial states are usually unstable. One can find in the 
Landau description of fluid turbulence the instability necessary for the 
excitation of turbulence. It is also known, however, that the excitation 
of turbulence in liquids is one of the most complicated problems. The 
same is true in plasmas. Many plasma instabilities are known; they can 
be divided into so-called kinetic and hydrodynamic types. For the first 
kind the growth-rate y is much less than the characteristic frequency of 
collective motions: 





Vkin 
oh) <i, (2.3) 


For the second kind the growth-rate and the frequency are of the same 
order and for so-called aperiodic instabilities the growth-rate is much 
larger. 

It is not obvious that the presence of instabilities leads to the excitation 
of turbulence, but it is generally the case. The important factor is the 
time needed for turbulence to develop. It is known from non-linear 
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optics (when large amplitude waves interact in solids, for example) that 
the first stage of wave mixing can be described as an instability of the 
generated waves. The instability is aperiodic, but this only describes the 
initial energy transfer. The full history is obtained from an exact treat- 
ment of three wave interactions. 

One can now ask if an aperiodic instability of a plasma leads to only 
one mode being excited. From the point of view of the uncertainly raised 
by the growth-rate y one should say yes, but one should really examine 
the stochastic behaviour of the exciting field. Let us introduce the time 
t needed for a plasma to become stochastic in its behaviour. If only one 
mode is excited it is obvious that the spread of energy to other modes 
is due only to non-linear interactions. The turbulent state is reached, by 
definition, when the energy is distributed over a large number of modes. 
The time needed to develop turbulence in the case of an aperiodic 
instability is, therefore, much larger than the growth-time of the instabil- 
ity. It is possible that the parameters which characterise the plasma state 
are changed sufficiently that the condition (2.3) is fulfilled in the last 
stages. This is the case, for example, if a low-density beam interacts with 
a plasma. In the first stage the instability of the beam is a dynamical 
one and aperiodic. The beam comes to a stochastic regime as the instabil- 
ity becomes kinetic. It is not known quite how this transition occurs: 
two possibilities are an explosive non-linear instability or the action of 
trapped particles. One can now not exclude the possibility that in the 
aperiodic type of instability there is a stage of its development when it 
cannot be considered a weak one. It is quite probable that the final state 
of stationary turbulence is weak, but it is not known if the statistical 
approach can be applied during this transition. If not, this means that 
although the energy in collective motions during the transition may be 
very high, this is not a strong turbulent regime. All of these questions 
could be solved experimentally with the present level of knowledge, but 
this has not yet been done. 

From the general statistical point of view, in an ergodic system the 
average over a statistical ensemble must be equal to the time average. 
This can also be checked. From the work of Fermi, Pasta and Ulam 
(1955), Kolmogorov (1955), and Arnold (1963) one can see that the non- 
linearity plays a significant role in the stochastisation processes. The 
smaller the difference between the characteristic frequencies of the system 
and the higher the amplitude of the field, the sooner the stochastic regime 
appears. In this sense, one could suppose that the best condition for 
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stochastisation is the presence of broad branches of frequencies of 
collective motions. This is, indeed, true for an infinite plasma as can 
be seen from Fig. 2, and is also true in finite plasmas, if there are modes 
with wavelengths much smaller than the size of the system. The presence 
of a broad spectrum means that the energy of the collective motions can 
be spread by non-linearity, at least in such a broad interval, and stochastic 
behaviour is quite probable. 


2.2. Description of Weak Turbulence by the Concept of “Elementary 
Excitations” 


Suppose that a weakly turbulent state exists in a plasma. It is usually 
said that weak turbulence has an advantage over strong turbulence in 
having a small parameter y/e (where y is the linear or non-linear growth- 
rate) in which everything could be expanded. The presence of a small 
parameter is a fact, but the conclusion that an expansion procedure is 
possible is mostly wrong! For example, strictly speaking one cannot 
expand the non-linear plasma equations in the energy of the turbulence, 
W, or the strength of the stochastic part of the field connected with such 
turbulence. We can give many examples of this. Only in an approximate 
way for some integrated value can such expansions be made, and there 
is a good physical reason for this statement as will be shown. 

If a statistical description is valid and the collective motions of the 
plasma are of a wave type, as in Fig. 2, one can consider waves as 
“elementary excitations”, generally called plasmons (Tsytovich, 1967, 
1968b). The interactions of plasmons can be treated by perturbation 
theory. For example, one can introduce the number of plasmons, N,, 
connected with the energy W, describing the turbulent spectrum by the 


relation 


_ @k)Ng 
W, = “ne 7 (2.4) 


where c(k) is one of the frequency branches shown in Fig. 2. (This differs 
from that used in quantum mechanics by a factor # (Dirac’s constant). 
This will not matter as we know the N, can be described classically and 


h does not appear.) 
One can introduce the probability of, for example, the decay process 


shown in Fig. 10, 
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kg 
Fic. 10. 
Decay of a plasmon (,) into two plasmons (k, and ks). 


Denoting the probability of this decay by w(ki, ko, ks), we can write 
down the balance equation that shows the growth of the N,, excitations 
due to the coalescence of N,,, N,,, and its decay into Ny,, Nx,: 


ON, 
“ar = | we ke, ks) [—Nx,(Ne,+1) (Neg+ 1) 
ak, dks 
+(Ni,+ Na. Nal aye 


dk. dk; 


~ (2x) . (2.5) 


= | w(k1, ke, ks) {Na Nae Nu Nag + Neue} 


The normalisation of the probability is to the phase volume dk ; d°k2/(27)8. 
Similar equations can be written for the processes shown diagrammatically 


in Figs. 11, 12 and 13. 


Fic. 11. 
Emission of a plasmon by a particle. 


(a) (b) 
Kia. 12. 
Two possible mechanisms of plasmon-particle scattering. 
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Here the solid line corrresponds to particles (electrons or ions). The 
process shown in Fig. 11 decribes the quasi-linear effect which can be 
called induced Cherenkov emission and absorption, first investigated in 
plasmas by Drummond and Pines (1962) and by Vedenov, Velikhov, and 
Sagdeev (1962). The process shown in Fig. 12 is an induced scattering 
which gives a non-linear energy transfer in the spectrum of turbulent 
waves. Sometimes it is wrongly called a non-linear Landau damping. It is 
not a damping at all, because it can be shown to conserve the number of 
quanta. This process was first mentioned by Pauli (as an induced one), 
but he did not consider the process shown in Fig. 12b, which is very 
essential in all turbulent plasma-energy transfers. In the case of high- 
frequency transverse waves the contribution of this effect is small com- 
pared to that of Fig. 12a, but it generally exists, if one considers the effects 
of many particles, even if they are only free electrons. For the ion-sound- 
wave interaction with the ions the process shown in Fig. 12 was first 
described by Kadomtsev and Petviashvili (1962), and the non-linear 
interaction of Langmuir waves was first considered by Sturrock (1957). 

The process shown in Fig. 13 is a four-wave decay process—a scattering 
of plasmons by plasmons. 


me OS 


Fic. 13. 
Two possible mechanisms of plasmon-plasmon scattering. 


After the balance equations are written down the probabilities can be 
found by the correspondence principle, that is, by finding the power 
emitted in the Ny, mode, if Ny, + 0, For the case of equation (2.5) this is 


ON«, a3 ky 
as Joa or dt (2x 


a3k, d*ky d®k 
= fv, ke, k3) Nu,Nx oka) ee 


an? (2.6) 
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One must then compare this with the work done classically, which ts 
Q=- [oe d°r. (2.7) 


This approach and the interpretation of non-linear interactions as 
induced scattering processes was developed in detail in papers by the 
author and collaborators (see Tsytovich, 1968b). 

It is possible to find all the required probabilities in the presence of an 
external magnetic field—for example, the linear and non-linear inter- 
actions of drift waves, and also of the modes shown in Fig. 2. This subject 
will be treated in full in the next section. 

Now we must return to the approximations needed to describe the 
turbulence by means of elementary excitations. First, one can easily see 
from the non-linear interactions described by Fig. 12 that the plasmons 
are driven in k-space and flow to smaller k (as was mentioned earlier). 
Such a process has a finite characteristic time, t,. It is apparent that 
the frequency of the plasmons will be uncertain with an uncertainty 
Aw ~ 1/t,. Therefore, there does not exist a one-to-one correspondence 
between w and & (energy and momentum of the plasmons) and for each 
k one can only have frequencies in some broad area around the dispersion 
curves shown in Fig. 2. Thus, the physical nature of the plasmons in 
the turbulent state is quite different from that in statistical equilibrrum 
because the frequency width is determined, not by the decay of the plas- 
mons, but by their non-linear interactions with other plasmons, which is 
dependent on the level of turbulent energy (which determines the number 
of other plasmons). 

Usually there exist some resonance conditions for the interactions 
shown in Figs. 11, 12, 13, which are nothing but the laws for conservation 
of energy and momentum during the interaction. For the process shown 
in Fig. 11 it is 

aw(k) = (ken), (2.8) 


where v is the velocity of the particle. 
For Figs. 12, 13 respectively: 
o(k) = w(ki)+o(k—k), (2.9) 
co(k1)—(kiev) = o(ke)— (kee), (2.10) 
co(k1)+o(ke) = w(ki)+o(kit ke — ky). (2.11) 
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Because the frequency w(k) has an uncertainty these laws only need to be 
satisfied within the accuracy corresponding to this uncertainty. Thus, 
resonance interactions such as (2.8) are not proportional to 


d[w(k)—(k-v)], (2.12) 


but to a broadened function of finite width. This point is essential, for 
example, for ion-sound interactions with electrons, when the broadening 
of the resonance (2.8) changes the interaction appreciably. As will be 
shown later, this effect comes not from the uncertainty in @, but is due 
to the effective turbulent collision frequency which broadens the reso- 
nance (2.8). The broadening of (2.9) is also essential for ion-sound inter- 
actions. 

One can see now that it is possible to derive the equations for plasmon 
interactions as an expansion in the magnitude of the turbulent energy 
only if one integrates over all w. This takes into account the whole area 
around the broadened dispersion curve of w(k) and neglects approxi- 
mately the uncertainty due to the turbulent broadening. Indeed, in such an 
approach it is possible to derive the balance equations given above from 
the Vlasov kinetic equation. Dividing the distribution function into two 
parts 


PaT OS Oa Se, (2.13) 


and similarly 
E = Evt4Estoch, —(Estoch) — 0, (2.14) 


one then expands in E%°**, constructing an equation for (Etech Byte 


which is then integrated over w. It can then be seen that the regular part 
of the distribution function describes the solid lines of Figs. 11 and 12 


and the J Wy, do = N,w(k)/(22)° describes the wavy lines. The process 


shown in Fig. 12b can be interpreted as a scattering of a plasmon by 
the oppositely charged cloud surrounding the initial charge. Thus, f"® 
does not describe real particles (they take part in both the free motion 
indicated by solid lines and the plasma oscillations indicated by wavy 
lines), but an elementary excitation, that is, it represents charges sur- 
rounded by polarisation clouds. 
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2.3. Problems of Correlation Broadening 


Suppose thata stationary state of weak turbulence is achieved; then the 
principal question is that of the structure of the correlation functions 
in such a state. To construct a theory of correlation functions it is obvious 
that an expansion in the energy turbulence, W, is not possible even ap- 
proximately. Indeed, if one tries to use an expansion similar to that used 
in deriving the balance equation one finds that, for example, the correla- 
tion function for longitudinal waves is proportional to 


1 | Sone 
e(k)  w—wa(k)* 
Here o(k) is the solution of the dispersion relation e(k, w(k)) = 0. 
This resonance does not matter for the balance equations because, after 


integration over w, only the imaginary part of (2.15) contributes. This 
can be approximated as 


(2.15) 


n8[co—o(k)]. (2.16) 


However, for the correlation function it is necessary to describe the 
whole correlation curve inside the resonance. AS was mentioned above, 
the width of this resonance depends on the energy of the turbulence. 
Formally, this divergence comes from using the Maxwell equation, 


e(nE, = i Soecteas i Foy Eg (#1 —H2) de, dt, 


n= {ko}, dx = dk do, (2.17) 


where the S, ,.,,. are the components of the non-linear current, to express 
the average (F,F2F3) in terms of the average of four fields (FE, F,F,£,). 
It is obvious now that to divide equation (2.11) by e(*) near the resonance 
is impossible. It is possible to construct a non-linear integral equation 
for (E, E,,E,,) by supposing that [w—«(k)]/o@ and W/nT are of the same 
order of magnitude. Thus, it is necessary to use a more precise equation 
then equation (2.17) including higher-order terms in the expansion in E,. 
A term proportional to W in ¢() arises and gives the turbulent renormali- 
sation of ¢(x), which takes into account the turbulent collisions. This 
kind of approach was derived in a paper by Makhankov and Tsyto- 
vich (1970) and in a review paper at the Bucharest Conference (Tsytovich, 
1969a). 
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The complicated non-linear equation for (E, E,,E,,) can be solved 
approximately if one supposes that in the terms containing the integral 
over w, 1/e(x) has no resonance behaviour and can be approximated by 
(2.16). In other words, it is possible to make the same approximation 
as in the balance equation, but only in the integrated terms. The final 
result given below in equation (2.18) for the correlation function has a 
force that can be compared with that found by the expansion procedure 
and given in equation (2.19). 


Lfe(x)+ 2" (%)] 





I, I, | S,, yy \? d(x — x1 — x2) d4 41 d' 2 
—= 2 PL OB h VO, SBN Se N OO hed Se te ee Tat 
= 32n [- fel te" (—w)] , (2.18) 
T[e(x) + e® (x)] 
[, I, | Sx, x1 sXe [? o(x— x1 — m2) dey d* 2 
a 2. 1°%2 <1 
= 32n fF Ga) » (2.19) 
where again 
x={k,wj, d'x = dk dw. (2.20) 


The difference is that the non-linear dielectric constant e™! occurs in 
the resonance denominator. The quantity «™!- is defined as follows: 


82i 


eh) =o [Banh db (2.21) 


82iS,., Hy Mm My Syn, x Hy 


(w—@1) e(x—41) G22) 





a aan ee 
where S and » are the components of the non-linear current: 
je = 0,.Byt | S12 Ea Eng 5(t 12) ey 
+ [2 segs Xp 9 Lv Losey Eng 5(H — 11 ~%2 — 3) U4) d4ug 4x3. (2.23) 


S and 2 can be found explicitly by the usual expansion of the Vlasov 
equation. It should be mentioned that the balance equations found from 
equation (2.18) using approximation (2.16) are practically the same as those 
found from equation (2.19), Equation (2.18) can be written in a form 
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which shows that J, is always positive: 


I, 


oe { Fry | Sx, nu nal? 5 = a2) bes Ulta 09 9) 


wo fe(x) + 2! (x)]? 


This necessary and reasonable condition is not always fulfilled by 
various approaches to the description of fluid turbulence (that is, strong 
turbulence). 

From equation (2.24) it is easy to see that the correlation curves have 
a Lorentzian structure near the resonance w = w(k), if one neglects the 
difference between and w and w(K) in all terms except ¢(2): 








Wik” 
Wie Lene 2.25) 
bo = Toa) +O] 
where 
nt. Im [e(e) +e" -(x0)] 
Ye ~~" Fee) Bem bcd (220) 
“by = ok) Ree 2.27 
@ (k) o(k) Be(x)/Oeo at : ( < ) 





Equation (2.27) describes the non-linear shift in the frequency. 

For weak turbulence yt! « w(k). If the turbulence is stationary, yt! 
may be small because of the non-linear compensation of the linear damp- 
ing and even smaller than the collision frequency. Therefore, the order 
of magnitude of yt" corresponds usually to the slowest decay process 
allowed by the conservation laws. As can be seen from the balance 
equations (2.5) the decay balance is not a differential one, but an integral 
one, the turbulent energy coming from different & values, ke, ks. 


2.4. Problems of Stationary Turbulent Spectra 


To find the stationary spectrum W, it is possible to avoid the detailed 
forms of the correlation functions and simply use the balance equations. 
Such equations were first used by Kadomtsev (1962) to find the spectrum 
of ion-sound turbulence in a partially ionised plasma. Then in some 
later papers there were attempts to use such spectra for a fully ionised 
plasma, which, as we shall see, is sometimes doubtful. The Langmuir 
turbulence spectrum was calculated by using the balance equations in the 
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papers by Pikel’ner and Tsytovich (1968) and by Liperovskii and Tsyto- 
vich (1969) and numerically by Makhankov et al. (1970). We shall have 
some special remarks to make later about both Langmuir and ion-sound 
turbulence spectra. is 

Here we mention that after finding the spectrum it is possible to calcu- 
late the correlation broadening. In particular, itis possible to use equation 
(2.24) to find the form of the correlation curve outside the resonance, 
that is, on the tail of the correlation curve. This is shown schematically 
in Fig. 14. 


Ik. 





Fic. 14. 


The correlation curve of Langmuir turbulence at low frequencies 
(do/w ~ ko/k). 


The whole of this curve can be found when the power input Q is known. 
The maximum in Fig. 14 corresponds to frequencies equal to the differ- 
ence between the frequencies of the turbulent waves. These correlations 
have no w, k-dependence. Their action on the particles gives the heating. 
Because the energy in the tail is proportional to W2 (see equation (2.24)) 
the next terms compared to the quasi-linear diffusion, which is propor- 
tional to W, in the expansion of the diffusion coefficient in £, are of 


TSY 3 
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the same order of magnitude. Both terms together describe the heating 
due to induced scattering. 

Now, if the stationary turbulent spectrum is known, it is possible to 
find the general form of (i) the rate of stochastic acceleration of fast 
particles, (ii) the scattering, amplification and fluctuations of electro- 
magnetic waves propagating in a turbulent plasma, (iii) the anomalous 
turbulent emissivity, and (iv) macroscopic characteristics such as the 
anomalous conductivity and anomalous diffusion. 


2.5. The Electromagnetic Properties of a Stationary Turbulent Plasma 


As was mentioned previously, for low frequencies, @ < v4? (oP is a 
characteristic turbulent frequency), the electromagnetic properties of a 
plasma for small signals or small perturbations are changed significantly. 
This means that the dielectric constant ¢ as a functional of W,, 


e = e(a, k, Wi), (2.28) 


cannot be expanded in W,,. One possible method of deriving (2.28) is 
to consider a small deviation of this stationary turbulent state due to a 
regular electric field perturbation, E", and expand all variables in the 
amplitude of E'°®, Collision integrals arise which describe the turbulent 
collisions. It is sometimes possible to expand the kernels of these collision 
integrals in terms of the turbulent energy, that is, to take into account 
the highest turbulent collision rate. Nothing is expanded in v/a, and 
in this way the essential change in « is found. Sometimes this expansion 
of the collision integrals does not work because a resonance denominator 
arises. Thus, a problem similar to that of the correlation broadening 
arises. The non-linear integral equation for summing the series in the 
turbulent energy for such collision integrals can be constructed. It can 
be solved and the collision integrals found. This give e(w, k, W,,) and 
the electromagnetic properties of a turbulent plasma can be investigated. 

The first approach to such problems was made by Rudakov and 
Vedenov (1964), who used a phenomenological approach with a Miller 
force, and by Gailitis (1966), who used the correspondence principle. The 
summation in the turbulent energy was done by Tsytovich (1968a, 1969b) 
and an investigation of drift waves in a turbulent plasma was given by 
Krivorutskii, Makhankov, and Tsytovich (1969). There is also an 
approach in which the turbulent collisions are taken into account 
phenomenologically. 


3. The Balance Equation for a Turbulent Plasma 


3.1. The Refractive Index for Waves 


The balance equations can be found in a simple way using the concept 
of induced processes. Let us consider a linear mode propagating with 
small damping, or small excitation, in an inhomogeneous plasma in an 
external magnetic field, so that 


o = otiye, (3.1) 


where yz < wz. AS we mentioned above, yg describes a linear or quasi- 
linear growth-rate or non-linear energy transfer. In the first approxima- 
tion, one can neglect all these effects if y < @, and say that equation (3.1) 
describes approximately a linear mode @ ~ w;. As we shall see later 
this wz is determined not by the exact distribution function f but by the 
part f£, which varies slowly in time. Thus, we consider w; as the frequency 
of anelementary excitation, and drop the superscript “reg”, but remember 
that f describes the distribution of the other elementary excitations—the 
“dressed” particles. Thus the dispersion law for the elementary wave 
excitations, the so-called plasmons, can be found from the linear disper- 
sion equation. One can also find the normal unit vector eg, which charac- 
terises different electric field components. It is useful then to introduce 
the dielectric constant ef which describes such a wave, or the refractive 
index n° = 1/(Re e’), where n° = k/wg and thus we have 


kK? = (wf) Re €°(wz, k). (3.2) 
Thus, wg satisfies the usual equation for transverse waves, but the polar- 
isation of the wave is in general arbitrary and is defined by ez. The 


question is how to find ¢’(@, k), if eg and the dielectric tensor ¢,(@, k) are 
known. The general linear dispersion relation is (c = 1) 


Y [k?5;;— kik; —w*ei(c, k)] (E,); = 0. (3.3) 


If the field has the direction ef one can put (E,); = (eg)£, and by multiply- 
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ing equation (3.3) by ej, and using 


(ef-e;) = 1, (3.4) 
we find 
[k? — we, K)JEX = 0, (3.5) 
where 
7 4 ; kez) (k-e?" 
e'(w, k) = (ef"); ico, k) (ex); + ea) f (3.6) 
One may say that to find €° is useless because to find ef one must solve the 
dispersion equation (3.3) for w(k) and put it in (3.3) to find the relations 
between the electric field components. However, e” describes the energy 
of plasmons and comes into all probabilities of non-linear interactions. 
Let us, for example, calculate the average energy density W of the field 
corresponding to the o-mode, supposing that the field has a stochastic 
nature. 


From the energy conservation theorem we have 
ow 1 oD 1 OB 
and putting 
f= ( Ee lt PR de, 





D; = | £1 (0, R) ER, (ef) a? deo, (3.8) 
B,. o i En, el _ lk A ei] Ex, os 
o ra) 
we find that 
ow 


oo = = a x | (i), Fi, k) wer); + (ek); 0’ elo’, k’) (ev); 
mid 


t 


Gu ti6)) eae (hee (eee) (et) ap wo ER, o! 
Www (Ww 


XMM TH dy Ae! (3.9) 


The result is symmetrised in x and x’. Integrating over ¢ and taking an 
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ensemble average, and using 





(Eko Ek, w) = [Eb (o+o') (k+K'), (3.10) 


one finds that 
|Ex, o [? 0 


gaa tan (€%(@, k)w?) dk do. (3.11) 





W) = 


In the approximation where the correlation broadening is neglected 
|Ex, wl” = aN (co — wf) +ayN_4 (wo +f) (3.12) 


(here wy > 0), we have 


— I ] é 250; 3 
(W) = a fo a Bo [w E (w, k)] | end k. (3.13) 
On the other hand, 
oN, d°k 
so that 
1 1 @ -1 
=-—|-7, -— {we (o, k ‘ 3.15 
01 5 [apt an Ms ay 


3.2. Quasi-linear Equations 


The simple balance equation, which describes the induced emission 
and absorption of plasmons by particles was first derived by Vedenov, 
Velikhov, and Sagdeev (1962) and by Drummond and Pines (1962). We 
shall write it in a general form applicable to a plasma in an external 
magnetic field and including possible inhomogeneities of the plasma 
perpendicular to the field. It is useful to start from the quantum de- 
scription of free-particle motion in a magnetic field and then consider the 
classical limit, or more precisely, the quasi-classical approximation. 

The particle motion in a magnetic field H, = H, is described by the 
energy spectrum of the Landau levels which in the general relativistic 
case are given by 


e = mct+cp2+c*p? , (3.16) 


30 THE THEORY OF PLASMA TURBULENCE 
where p, is the particle momentum component along the field and p’, is 
pi, = 2nh |eHg|, (3.17) 


where n is an integral quantum number characterising the energy levels 
and # = h/2zx, where fA is Planck’s constant; 


rt 0 bey ae (3.18) 


{n the quasi-classical limit the spin does not change and can, therefore, 
be neglected. The quasi-classical limit corresponds to a very high n such 
that n + oo, — 0, while nh is kept finite and p, is the particle momentum 
component perpendicular to the magnetic field. 

It is also useful to choose the gauge of the magnetic vector potential 
as in the first Landau paper: 


A, = Hoy, Az= Ay =9, (3.19) 
Hy = [TAA]. (3.20) 


The momentum component p, is conserved and the energy spectrum 
is independent of p,. In the quasi-classical limit p, determines the y-coor- 
dinate of the centre of the Larmor orbit: 


¢c 
y ae? Dx. (3.21) 


The probability of emission of a o-plasmon by a particle gyrating in a 
magnetic field (or in the quantum case by a particle in one of the Landau 
levels) is denoted by 


w(k, Pz Pr, v). (3.22) 
This probability depends, as written, on the momentum of the emitted 
plasmon &, on p,, the initial energy of the particle (or more precisely x, 
the integer characterising the Landau level), and on the final n’ in which 
one finds the particle after the emission. Instead of the variables n and n’ 
in equation (3.22), we use p, and », where 
y=n'—n, (3.23) 
and » varies in the range 


ya—oo,..., —1,0,1,..., te. (3.24) 


THE BALANCE EQUATION 31 


The distribution function f of the particles depends on p,, p,, and n: 


f =f (Pe Px N). (3.25) 
Expressions (3.22) and (3.25) are supposed to be normalised to the 
phase volumes dk /(2)? and d3p/(2x)'. 
The balance equation for the waves can be written in the form 


dN d 4 o o 
ae = ‘3 { o: [(Ne + 1) f.(p:, Px, n)—Nx f.(p.—hk,, px—hk,, Aa— ”)) 


xwe(k, Pz Ps »), (3.26) 


where a = e, i corresponds to electrons and ions, respectively. The term 
with the minus sign occurs because of detailed balance arguments and 
the conservation of momentum in the emission process gives 


P: = Pr—hk,, px = Px—hkx, (3.27) 


where the prime corresponds to the state after emission. In the quasi- 
classical limit Ak, « p,, hk, << p,, and »<n. Expanding in k,, k, and 
y it is useful to introduce classical variables for the argument of f: p, 
instead of n and y instead of p,. Thus, instead of k, Of/Op, we write 
—(k,c/|eHo|) (Of/Oy), and instead of » Of/On we write Ap’, Of/Op’,. From 
(3.18) one finds that 


Ap’, = 2vh|eHo| = hwyve. (3.28) 


Thus, equation (3.26) can be written as 


aNe i . 4, f(Pes Pas Y) 
pe ee d (2xy8 Nyw%k, Pz, Ps) Mj aaa) Raa 
+ wk, Po Ps”) FAP Pi, | (3.29) 


The first term describes the induced process, the second the spontane- 
ous emission, and 


Ai = Pz> Ac = Pi, As =y), (3.30) 
— OHa? 


c 
MA, = kz, Ake = , Ads =-k,——, 3.31 
ah : [effo| Ga) 


v= Es, (3.32) 
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Similarly, one finds the balance equations for particles: 


A Pz P1>)) 
oo de. = 35), 0% "6i; Oe epics N+D py Ai fl Pz Dis Y)s (3. 33) 


where 


oo 


- § = aa AA; AAW, Pz, Pi» YNk , (3.34) 


+00 
4; = Di | (22 Aiwy(k, Pz Ps v). (3.35) 


The first term of (3.33) describes the induced process and the second 
one the spontaneous emission. The d/dt operator on the left-hand sides 
of equations (3.33) and (3.26) describes the free transitions of particles 
and plasmons: 


o -7 +(o-v)f, (3.36) 
t 
dN IN, 
a = By t ers VINE» (3.37) 


where v,, = dwy/dk is the plasmon group velocity. In the case of an 
inhomogeneous plasma, equation (3.37) is sometimes rather more com- 
plicated and contains an additional term 


— (ve: a) (3.38) 


3.3. The Probabilities for Plasmon Emission 


Let us consider the power Q emitted by all the particles in the limit as 
Ng —~ 0, when only the spontaneous emission is essential in equation (3.29): 


o AN, d®k 


QO = | Gn) = oe O,, (3.39) 
i ‘dt (0) \yo.so J (20) Pz, Pi» WEp> ; 





where 


he Cr d®k 
= Y OW (k, Pz P11, v) (2m) (3.40) 


y= —00 
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is the power emitted by a single particle. To find Q, in the classical limit, 
we shall use a procedure similar to that used by Landau for calculating 
Cherenkov losses, that is, we will calculate the work done by the field 
E on the current j: 


C= { (j(r, t)-E(r, t)) dr. (3.41) 


The current produced by a particle spiralling in a magnetic field, 
j = ev(t) 6(r—r(2)), (3.42) 
where o(t) and r(r) are simple periodic (sine or cosine) functions of time, 


has Fourier components 


dies >is aa x Tye)" Se —k,v,—VOHa); (3.43) 


ro = {Xo, Yo, Zo}, Xo and yo are thus the coordinates of the centre of the 
Larmor orbit, z) and ¢ are the initial z and phase ¢ of the particles at 
t = 0. The components of the vectors I" in the fraine k, = 0 are 


(i). = v1 15; (Za) 5 (Li)y =v Si (@a)s 
(TD, =od(e); = ths 9, HPA aay 
WHa E 


The J, are the Bessel functions of order v. The field produced by the 
current (3.42) is found from the Maxwell equations (see equation (3.3)) 


[k26;; —kikj—we;,(w, k)] (Ex) = 4tion( j,i. (3.45) 
Because we are interested in the emission of the o-wave, we must set 
(E.); = (ek) Ex (3.46) 


in equation (3.45). Multiplying (3.45) by e%* and using the definition (3.6), 
we have 


ES = Asien at is Ho. (3.47) 


Expanding j and E in equation (3.41) in Fourier series and using 
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equations (3.46) and (3.47) we get 
ak dw do’ ° ny o oF (o-oo 
Cia CE FE otge Stil ie ex) (ef aj, Je °™ 
— wre 

4ni pak —i(a—w') +i(v—v') yo * oo o* 
7 (2m) Y, Pes } ° . i (Finel) ee 
x per agg 8O—Kevs— rps) 8(0'—kev,— Voy.) ak do doo’. (3.48) 
We substitute here the exact expression (3.43) for the current j,. 


Averaging over the initial phase qo, 


1 
cae i(v—’)}@o = ; 
an | e! doo = 4,,¥, (3.49) 
and noting that only the imaginary part of 1/(k?—qe°) contributes to 
equation (3.48) and that 
Poot © ni 5(k® —we’) 
nj WO og) + 5(w +04) 
GIO) (wer) 





, (3.50) 





w=o8 


one finds an expression of the form of equation (3.40) and can immedi- 
ately find the required probability 
2 | (Meez)|” Swe —k0,— Onn 
Wilh, Pes Pa, ») = (2m) — oie other te) (3.51) 


w= Ig 


Using this probability in equations (3.33) and (3.29) gives us general 
quasilinear equations valid, for example, for a relativistic plasma. The 
equations describe all processes of interest, including transverse waves, 
synchrotron emission, synchrotron instability, spontaneous processes, 
and they take account of plasma inhomogeneities. 

In the case of longitudinal oscillations, it is useful to introduce the 
longitudinal dielectric constant in a way slightly different from equation 
(3.6): 


LE ao = 3 Ls (3.52) 


E 
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The dispersion relation (3.5) then has the form 
el =0. (3.53) 
Using the fact that e, = k/k and that 
0 
2% a 1ye 7 
wey a = (a,) op 


we can write down a simple expression for the probability for the emission 
of longitudinal plasmons: 


“ (3.54) 


= Wp, 


deo 





w(K, Pes Dis ¥) = (2x)? — al eee a (3.55) 


wo=w! 
k 





3.4. The Case of Non-magnetic Particles 


Even with a magnetic field present, in some cases the Larmor radii 
of the gyrating particles are much larger than the wavelength, that is, 
2, >> 1.7 In this case, a large number of cyclotron harmonics are involved 
in the interaction and the probabilities which we gave above are useless. 
In this case, one can consider the opposite limit, by supposing that 
during the emission the particle does not gyrate, but moves rectilinearly. 
One can use as the particle current the expression 


= Ge = hae r) 8(o—(ked)). (3.56) 


The quasi-linear equations have the same form as in the absence of a 
magnetic field: 


af a 
oe ane "Se Le Aif, (3.57) 
Diy = | cae w(k, p) Nk; ik; > (3.58) 
dk 
Ai = x8 kw<k, p), (3.59) 
aNg d°p 


E a 
nN) ape (te 3p) +) Gems 090 


t We call such particles ‘non-magnetic particles”, 
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The probability can easily be found from equations (3.40) and (3.56): 











* (ok —(ke0)) (3.61) 


w(k, p) = (an) © - ee 





o 
D= Oy 


Note that equations (3.57) to (3.60) are also valid for a relativistic 
plasma or for relativistic particles in a cold plasma. 


3.5. Examples of the Quasi-linear Equations 


In the case of the Rosenbluth—Post loss-cone instability (Rosenbluth 
and Post, 1965) the ions are non-magnetic and one can use equations 
(3.57) and (3.60). Since k, >> k,, we have 





1 dNk 2 
AN« _ = Sawaal 8(col, — (kev) 9 (Kivi) Of d°p, dpz. 
Ne dt k*(de8/ dw) VL Op, 
oes (3.62) 
Introducing 
fle.) = gpate | £° be (3.63) 


one immediately finds Galeev’s quasi-linear equations for this instability 
(Galeev, 1967): 


josey 


1 dNi POY) " fles) dey (€1) de, 
yl ~ 24/(2) RK (Be,/80) 3/2 (3.64) 
Ny dt — 24/(2) kk, (de,/0m) J. et 
where 
I 1y2 
El,o= ini E1,0 << E11 min> (3.65) 
and 
s if E— <= E) min» 
3.66 
wale. y SE.) if @ > e) min» woe 
Ofer) 8 ( 1 Ofer) 
ee ee, (a: ‘ a Ce) 
(coy)’e °Ne AVAL 
— 3 eh AE 
a ‘f GT) PAK, (OHJ00)) oy 
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Another example: if the Larmor radii are very small (z « 1) one can 
use for the probability for the emission of longitudinal plasmons the 
following expression, which is easily found from equation (3.55): 


(2n)°e? O», 0 


Wk, Pes Bis) = Spar igay)  Moke—ker:). (3.69) 





| 
o=m 
k 


This immediately gives the quasi-linear equations for drift waves, first 
found by Galeev and Rudakov (1963): 


dN, = e2N, ky 
“dt ~ xmk(delJOm)| _ <etea mee 35 Olek ker) 


wre 











(3.70) 
A =| an(e. 2 ae ky -) (a 6(coh—k kw, Nx 
kz Ov: One Oy) ak?m2(Ge!/00) | 
lo =o} 
0 k, 0 
(eden a) oan 


It is evident that the method we have described can be used to obtain 
rapidly results which would normally require a substantial amount of 
work. 


3.6. Non-linear Plasmon-Plasmon Interactions 


Let us consider the simplest non-linear wave decay process: 
o>o'+o". (3.72) 
We shall denote the probability of this process by 
we" (k, ki, ke). (3.73) 


It is easy to obtain a balance equation similar to that of Section 3.2: 


dNy 


d*k, d® ke 
“dt , 


Day (3.74) 





= i we” (k, ky ke) {Ni Ne, — NENG. —NENz} 
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To find the probability let us again use the correspondence principle. 
The emissivity of o-waves in the limit as Ng > 0 is 
ov = w;, aN ak 
~ dt np 


dk dk, Ake 


= [ owe (ky kay ka) Ni, Nig Gg (3.75) 


This emissivity is equal to 
QO? =—((jrE)), (3.76) 


where j is the non-linear current excited in the plasma by the o’- and 
o’-waves, and E is the field produced by the current. 

We recall the general expression for the second-order non-linear 
current: 


(ii = ps { Silo, #1, Ho) 8 41 — 20) Ej, x, El, xp T41 d4x2. (3.77) 
Js 


The S,, can easily be found from the Vlasov equation by means of an 
expansion in powers of E. In the present case equation (3.77) becomes 


(ji = 2 f Si(x, x1, 42) Eg, Eg, 6(¢—x — x2) d*x; d*x2. (3.78) 


The fields E can be found from equation (3.47). Substituting for j into 
equation (3.76), we obtain 


a abi 2 o* onl. 
) C(egejnr”)) (ek ein’) 





O =—-4azi 


exp oe ear (3.79) 


We now use equation (3.50) and (3.78) and average the result by using 
equations (3.10), (3.12) and (3.15). Comparing these results with equation 
(3.75) we see that 


wi?" (eh, ka) = 32o(2n)" 8(k— hi ka) 8(o% — eof, — 0%) 











> (w4,)” (we, )? | Scor(ok, k, wy, k1, wi,, k))? (3.80) 
bus o Ow Orem, o’ OW2 O2E se ae 
wry lo =op = OK 
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39 
i,j,l 
and 


(3.81) 


(3.82) 
From this general formula it is easy to find many of the results of other 
authors. 


3.7. Non-linear Plasmon—Particle Interactions 


To a first approximation such an interaction is due to the induced 
scattering of plasmons by particles 


ota zo'+a’. (3.83) 
Let us denote the probability of this event by 


we (k, k’, Pes Pi, »)- (3.84) 
The balance equation is very similar to that for emission except that 
the momentum change is now 


P, = P,+k,—k,, Dy = Pxtki—k, (3.85) 
We have 
dN dp, a°k' 
= a Y [ae (Qn si We (ks k’ » Dz Px v) 
x [Ne (Ne + 1) fp: Px, n) 
— (Ne +1) Nef pet ke — key Put kx— ke, n—v)) (3.86) 
In the quasi-classical limit we obtain 


te d'p d*k P f 
= Pe Se (On NENG WE (kK, Pas Pas) ARS pg felPer PY) 
t= 0 d3p dk 
+X, 2% | Gaye 
“=f d8p! dk’ 
a,0° ¥=—o0o (2z)' aa 


Nyew, (k, k’, Pz Pi> »”) f.( Pz, Pi» y) 


— New? (k, k’, Pz Pi> ») (pz Pi, y), 


(3.87) 
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ky—ky 
Vy =k,-k:, Ady = --2--- , Ags=- se. 
: an 3 Teel 





(3.88) 
In the limit as N, > 0 equation (3.87) describes spontaneous scattering: 
Tt, Dat NEWER RS Pos Pi Mh B89) 


This can be used to find the emissivity, as follows: 


aA d3p aa’ 
=| sso". (3.90) 
aa oe ak 3k a , o 
Qn = _ (2xr)6 Nawo(k, kK’, Py Ps VOR (3.91) 


In order to find the probability it is necessary to know the current due 
to emission. This current has two components, corresponding to the two 
possible mechanisms of scattering, that is, the usual Compton scattering 
and the non-linear scattering (see Fig. 15). 


ra (a) (b) 


Fic. 15. 
Two possible mechanisms of plasmon-particle scattering. 


Often these currents cancel one another as, for example, in the case of 
the scattering of Langmuir waves by electrons, or ion-sound waves by 
ions, and in many other circumstances, 

The current due to Compton scattering is easily found. It is necessary 
merely to compute the small oscillations of a particle in the E°-field 
(c = 1): 


where 
He = | RAED 3.93) 
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In first approximation the particle moves either with constant velocity 
or spiral in the zero-order constant magnetic field. If this trajectory is 
substituted on the right-hand side of equation (3.92) the small oscillations 
of the particle in the field can be determined. The current due to these 
oscillations is easily found from the relation 


j = ev(t) 6(r—r (0). (3.94) 

Although in the relativistic case the calculations are sometimes tedious, 
the general result is usually quite simple. We write it in the general form 
. d‘x' 
bm dd. | Oy 


a,0° v=—oo 


Ag, x, 9B EOF" 8 (cy —00' — (kz —k 02 — POH), 


(3.95) 
without specifying the coefficients 7. 

The current due to the non-linear scattering is also easy to find. 
Indeed, such a current stems from the simultaneous presence of the EY- 
field and the field produced by the particle, neglecting the small oscillations 
of the particle in the E?-field. 

From equation (3.77) we thus obtain: 


Gi=y 2 | Sih, x, #—X') Ej Ew 1 dx , (3.96) 
a,o Jj,l 


where 
» [k?6,;— kjk; —we,(o, k)) (ED) = 42io( jn); , (3.97) 


and (j2); is given by equation (3.43). If one introduces the inverse Max- 
wellian operator by 


py TT, (x) (k?6i; —kik; —we;,(x)] = 45;, (3.98) 
one finds 
Ee, i= r 4xi(w — —o@ ’) 2 IT .(%— HX y, ba Joa Fhe x, s 


exp Wk “yy *ro) + a 6(w —w' —(k,—k3)v,— POH) 3 
(3.99) 


I” differs from I’ (see equation (3.44)) by a rotation of the xy-axes such 
that k, # 0. Thus we find 
+00 
eee [a 1 Beh 1+! 5(cy — cn! — (kp — ks) — vO H2) dx’, 
(3.100) 


v= — 00 
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(an: i y 82tie(w —w') 


Dey (eds Sule ', x2) Ta ea Ta 
SIiAS 


(3.101) 


The emissivity of such a current can be calculated in a manner analo- 
gous to that in the case of plasmon emission, Comparing the result with 
equation (3.91) we obtain 


Cn Alok bok, HIE 


Bey 8) 





w9(k, k'. Pes Pi. ¥) = , (3.102) 








° @ 
Ws oe o 0 ener 
where 


Ago = ¥ (e4): (Ai + AP"). (3.103) 


It should be noted that the inverse Maxwellian operator represents the 
virtual wave in the non-linear scattering diagram. Often such a wave 
can be considered to be longitudinal, in which case 


1 kik; 


we! k? ° 


IT, = (3.104) 


where ¢! is defined by equation (3.52). If, in addition, the non-linear 


scattering dominates (as it does for the scattering of high-frequency 
plasmons by ions) a simple expression for the probability results: 


6 2) olf 1 iy rot V ry) 2 
wi (k, k', Pr Pis v) ey 32(27) e | S'(ox, k, Op’, k » Ok — Ok, k—k )| 





= oo alk—k' |’ |e\(w, —oy, k—k’))? 
ei (Za) (cox — we —(kz —k;)v,— VOH) , (3.105) 
(wt? 2! i | 
i Ow 7 inp wt Ow’ sd a =o, 
where 
er ({A —k'], +01) 1 k; iki jkay 
Zag = sce an S' = =e kky ks Sij- (3.106) 


An equation for the change in the particle distribution function 
resulting from the scattering can also be written down, It has the same 
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form as equation (3.33) with 


+00 d*k d*k' 


Dy = Xd agp 14 Adiweke R, pes Dis NENG, (3.107) 
te 0 dk d®k' 7, ay 0) alot , 
A, = De eis - “Qn (Ne — Nz) Aliwe(k, k >» Pz Pi» v). (3.108) 


If the particles are non-magnetic (or if Ho = 0) then A; = p,; and 
AA; = k;. In the circumstances in which both the interacting and virtual 
waves are longitudinal, the probability is obtained from (3.105) by making 
the substitutions: 


PR 1, (3.109) 
(og — wh — (kz — k:)v2— vO) > (oe — OH —([k—k']-0)). (3.110) 
For example, in the simplest case, when Langmuir waves are scattered 


by ions, the probability is 


(2a)? Meek)? (ch — ol — (Uk —K']-0)) 


y Vaso SE A 
wi(p, k, k’) = new?, Kk? 1+T,/T; 


(3.111) 


The method described above enables one to obtain and use the probability 
for the process of interest. 


4ae 


4. Turbulent Collisions and Resonance Broadening 


WE EMPHASISED in the earlier chapters that the balance equations are 
only an approximation and that turbulent broadening of resonance 
interactions plays an essential role in the turbulent state. One can say 
that the broadening of resonances is due to turbulent collisions. The 
derivation which we have given of the balance equation shows that the 
particles emit and scatter waves in a linearly stable or unstable plasma 
as if there were a high level of dissipative processes. The point is that we 
used 


1 = 2 ay2eo 
Im wre = 1(k Ore ), (4.1) 


which is only true in the region of transparency where Im (e°) is small and 
positive. This implies the use of retarded potentials in calculating the 
waves emitted by particles. In an unstable plasma when the linear part 
of Im (¢’) is negative, one must use equation (4.1), but not the expression 
with the opposite sign (advanced potentials), because the linear part of 
Im (e’) has then no meaning for plasmons. It was already emphasised that 
the plasmons are driven in k-space by non-linear interactions which play 
an effective role as collisions. This means that the denominator in equa- 
tion (4.1) depends on the turbulent energy which can, therefore, not be 
used aS an expansion parameter. Only in the balance equations where 
the broadening due to the turbulent collisions is neglected can one use 
equation (4.1). Also, the transition to the turbulent state involves non- 
linear interactions and thus a more precise description of the turbulent 
state is needed to include turbulent collisions. 

Maybe one can say that for the balance equations it is an academic 
question, because it is obvious that in a statistical description one can 
use a transition probability approach and that the increase in entropy 
needs a retarded potential. For example, it is obvious also that binary 
collisions in a plasma could be described by collision probabilities or 
cross-sections and the exact derivation of the collision integral from the 
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Bogolyubov scheme can be regarded as an important, but academic 
problem. However, this is not so because Bogolyubov’s scheme is the 
way to prove the validity of the physical picture used. 

Let us emphasise then, that the question of turbulent collisions is not 
academic, nor is the treatment of the Landau poles in the balance 
equation, but that they are essential practical questions for ion sound, 
for example. The study of this problem by Rudakov and Tsytovich (1971) 
shows that the induced scattering of ion-sound waves by electrons be- 
comes negligible compared to the quasi-linear collisions. Therefore, not 
all possible diagrams should be taken into account, even in the balance 
equations, and a full statement of which diagrams are the most important 
will be given later. Let us only remark that the lifetime of the elementary 
excitations depends on the turbulence energy and that they are physically 
quite different from excitations near statistical equilibrium. 


4.1. The Balance Equation found by Statistical Averaging 


Let us first find the balance equations exactly from the Vlasov and 
Maxwell equations using an expansion in the stochastic field amplitude. 
For simplicity let us consider waves in an unmagnetised plasma which 
are not resonant with the particles, that is w ~ (k-v), and start from 
the equations 


Be 1 +o: Ufa )+e(E- a) > 0, (4.2) 
dp 
(9-B) = 4n Ye, | fens (4.3) 
We divide the distribution function and the electric field into regular 


and stochastic parts as follows: 


fe = fie faroch, — (fstochy — 0, (4.4) 
E = Es + Estoch | (Estoch) = (0. (4.5) 


Thus, by averaging equations (4.2) and (4.3) and substracting from the 
exact equation the averaged one, it is possible to find two sets of coupled 
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equations for the regular and the stochastic variables (f = f'°?): 


Fa townte(em-Z) = ~e( (Bw (4.6) 


Of stoch of opus 
yee — + (orate) +e( Ester os te(Eme ) 


ot 
tel (pone 0) e ( (Be we ie. Ga} 


3, 
(7- Ets) = 4n Ye, | te ‘oe (4.8) 
(V+ Estoch) = 4 Y c. | pte one (4.9) 


Let us consider for simplicity the case E™® = 0, and expand in Esto, 
In the Fourier representation equation (4.7) splits into a series of equa- 
tions: 


ifo—(k-v)] fuer ( Bae. ») =0, (4.10) 


—ifo—(k-v)] ft" +e | | (Ex , te) 


2 eaten 
-((ee*-! Eee “) | dees) a dx, = 0, (4.11) 


— if — (kev) FOP +e | ( Estoch. oft. te) 


Op 
cl 
2 ( (Ee —)} B(x — 21-09) dy d4g = 0, (4.12) 


Decne ee eer een e eee e renee eee renee ere e eee e sees essere ee ssnr esses eer eee nas eeeennees 


Here the /" are for simplicity taken to be time and space independent. 
The generalisation is simple but this approximation is appropriate when 
Jf*® varies slowly in space and time. 

The first step is to substitute from equation (4.10) for f°", given by 


De e Estoch of 
PP = (*-35)} (4.13) 
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into the “collision integral” of equation (4.6) to get the quasi-linear 
equation: 


h 
= “I (Ba( Estoch (z oe -)) ei(ths +kho}-r)— ilar +eoo)t d'x, Ax 


0 bee ell{ky +h2}- ry- foi tony of 
a ( Een I sega 
ra oe kake eee op 
0 of 
= Lar (Dia) 4.14 
d api ( 4 ap) nn 
where 
kiki; 
Dy = me | a | Eas |? (1 — (hs -0)). (4.15) 
1 
Here we have used the usual expression for the average of the stochastic 
longitudinal fields, a minus sign arising as ky = —ko, 
(Ei Eee) = —| Ege |” 524+ 2). (4.16) 


The imaginary part of [w—(k-v)+i6]"1 has been put equal to 
—n6(w—(k -v)). The essential question is how did we decide on this rule 
for dealing with the singularities [o—(A-v)]~!? One could, of course, 
introduce a small collision term in the initial equation (4.2), but this is 
wrong because the turbulent collisions are much larger than the binary 
collisions. From expression (4.13) one can see there is a divergence near 
@ = (k-v) and one concludes that near the resonance one cannot use an 
expansion in Es'°** and the higher-order terms in E3'°* in equation (4.10) 
are essential. As we shall see this gives a broadening of 6(w—(k-v)) 
in the quasi-linear equation (4.15). Because this function has an appre- 
ciable maximum near @ = (kK -v), one can use equation (4.15) as an approxi- 
mate expression for the diffusion coefficient. This can only be done for 
integrated functions such as the diffusion coefficient, but not for the 
stochastic part of the distribution function f%'°, so that equation (4.13) 
is wrong near the resonance. 

Let us leave this problem for the moment and consider the case where 
@ # (k-v), when expression (4.15) is approximately zero and equation 
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(4.13) has a finite source term. An example of this is found in Langmuir 
turbulence, when the phase velocity of the waves is much larger than the 


thermal velocity. Substituting expression (4.13) into equations (4.11) and 
(4.12) one finds the non-linear charge density 


(2) Stoch (2) ap 
aeas } fa Ons 


are 2! eee (x 2) 5 ese (x 2) 
ky ko(2x¥ [o—(keo)] \' Op) w2—(ka-v) \" Op 
x [ Bstochesoch hs (Esch gsioch) &(x— x1 — x2) dx, ds» 


“ of he Pa 6 Sia eee ( Estoch pstoch] (x — x1 — x2) d’xy d's , 
(4.17) 


a3 
_ stoch(3)_@°P 
=e mel fe OF 


stoch stoch F ae stoch / ,Stoch y»stoch 
ae mir masta Ein Eg — En Eg Es 


— (Enoch pork p08] 8 (4¢— 2) tg) d'xe1 d'xy d'xz, (4.18) 


(3) it | 2 I (x 2) See eee 
Gomrnne =H J Gas Todo) Mt ap) @—ai— k= hilo) 


x (He: wy) saa (4s: x) Bo (4.19) 


One can now write the non-linear Maxwell equation (4.9) in the form 


ik 6, FeO — ax | 0... wl Es nai | ea = ( a Bstoch) 
x 6(x—m1 x2) d*x1 d’xe 
+4n fe Oe atu [ ee | ipa jaa Ee if ¢ Deieaay porn 
_ ( | alas 1 elas pen 
X 8% m1 — Hem) doer dhieg dies, (4.20) 
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where 

(i v=) 

_ 4ne? " @p _@p 

& = +s wo (Kev) 415 (2x) 7 (4.21) 


Note that «, contains the average regular distribution f"**, but not the 
linear stochastic distribution f%'°", so that equation (4.21) does not give 
a linear dielectric constant. Also it is necessary to mention that expression 
(4.21) was derived from equation (4.17) which has no exact meaning. 
But expression (4.21) is an integrated value, and use of Im [w—(k-v)+i6]7! 
= 76(w—(k-v)) that leads in equation (4.21) to the quasilinear 
growth-rate can be shown to be appropriate. Now w 4 (k-v) and the 
imaginary part of ¢, is zero. Using equations (4.20) and (4.21), and 
averaging over an ensemble, one has 


9 - t , 
by | | emacs |? = Dues: 3 ( sagas Esch Bo) (x — x1 —x2) da at d'x, d'xs 


+ 4ni (3) [ ( pon | eae | Saige Boch 


k Ox, My, Me, XQ! 
= ( ‘Estoch Eo) 4 Estoch Eso) 
x 5(x —#1—x#2—- 3) diy! dx, axe dx . (4.22) 


This equation is written down for the case of stationary turbulence. If the 
turbulent state depends weakly on time and space it is useful to define 
| Estoch |? by the relation 


Fa = | @ croaaesd raga elk tk 1 oto db ; (4.23) 


which coincides with that used above, if the turbulence is stationary. 
Then the equations for (¢,—¢,) (E°PE°) lead to the term 


which occurs in the balance equations. Let us emphasise here that the 
weak dependence on ¢ and r in the treatment given now is not a result 
of the smallness of the quasi-linear growth-rate, because the non-linear 
interactions may well lead to a stationary state. 
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We shall continue to consider the case of stationary turbulence and 
use equation (4.22). If the turbulence is weak, the correlations of the 
fields are rather weak so that one can put in equation (4.22) 
< | as Oia | oie eae a / Cpe Bae) ¢ | haa Datta 

rs ( y otaaaa Conia < Sib Ey a ( peo Ee) < Boot Eso : ( 4.25) 


The last term of equation (4.22) is, therefore, 
Ani a) [o®,,, iors r oo er | [Ene 2 eee r d' x, : (4.26) 


To find the average of the three fields one needs to express it in terms of 
four fields. This is possible by using the first term of equation (4.20) 


toch Ani 2 toch pstoch toch pstoch 
ES och _ a | 0? Vda ” [ES oc ES ‘oc et eee sa y] 
od 


X b(4—H1—H2) d41 d4x2 . (4.27) 


Let us mention that this expression has no meaning if e, is close to zero, 
as can be seen by analogy with our discussion of equation (4.13). But if 
it is used, one finds the first term of equation (4.22) in the following 
form: 


4 Gar stoch ;2 a ~(2) |ES 
|ES LE re Op ante Se soea Mage Ne i eee 





4 4 oy (2) ~(2) | more ee P 


kee_ On, x1, x20 —x, —m1, — x2 b(x Hy —xe2) dn, ‘xe ’ 


(4.28) 


where 6 is the non-linear charge density coefficient symmetrised in x; 
and xe: 


BP vce = FOP vse + Orne, o)* (4.29) 


One may assume that the first term of expression (4.28) will cause no 
trouble, because «,*,, occurs under the integral sign, but the second 
term goes to infinity when «, - 0. This inconsistency comes from equa- 
tion (4.27). However, let us continue the considerations by neglecting 
the term (4.28). This is possible, if the turbulent oscillations cannot even 


roughly satisfy the decay conditions required by the 6-function in the 
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second term in expression (4.28). This is exactly the case for Langmuir 
turbulence because 2,, # Wp, and the frequency of the waves lies close 
tO Wye. Now equation (4.22) can be written in the form 





(e, ugh nl. ‘) |Esoce — (4.30) 
where 
a = i Em Earn? doe, (4.31) 
_ 4ni 7 @) (3) (4) 
2x = ~ zo (ek many mb Oxy my — oy +4 Rk 
om a ‘ 
a Paget eae (4.32) 


The quantity 2® contains e® and the quantity 2°) similarly op. The 
term —Im e2"- /(0e/O)|,—., = yo" can be called a non-linear growth- 
rate and indeed describes the energy transfer of the turbulent waves. 
Because w # (k-v) the imaginary part of the’ first term of 2, ,, can 
easily be found from equation (4.19): 


Im 2, = as — 


~ (Hp) = ce (* “Wh 
int | Too (ean ) 


X 6(@ —@, — ([A — 3] 0) - (4.33) 


a aa 

By expressing |E,,|? in terms of N,, and calculating yh we easily 
find that equation (4.33) is exactly the same as the one we find for the usual 
Compton scattering of particles. If one examines ne next enn ae). 
and calculates the total imaginary part that exists in ¢, or o) Ma, x—ny and 
OP ice -x» one finds the non-linear scattering and the interference of 
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the non-linear and the Compton scattering, so that—as must be the case— 
the scattering cross-section is |M1+M2|? where M, and Mz are, re- 
spectively, the matrix elements of the two scattering processes. Thus the 
result is just the same as found in the scheme described in the previous 
section. 

To finish this kind of consideration, we return to the problem of stochastic 
heating. Supposing that for the turbulent oscillations w 4 (k-v) for the 
thermal particles so that the diffusion (4.15) seems to be absent. This is 
not quite true because the non-fulfilment of the condition a, = (k-v) 
in the centre of correlation curve does not mean that this condition 
cannot be fulfilled far away from this centre, and this indeed is possible. 
Jt is interesting that it is possible to find the form of this correlation 
tail by using the secular term (the second term in expression (4.28), which 
we have neglected in equation (4.30) or near the centre of the correlation 
curve). Indeed, far from the centre ¢_, has no singularities and the 
terms that are proportional to |£,|? become negligible because |£, |? is 
small. If one takes into account only the second term of expression (4.28) 
and substitutes it in the diffusion coefficient one finds that 6(@—(k-v)) 
is replaced by 6(w@—«@,—([k—k,]-v)), and the whole term is proportional 
to |E,, |? |E,, |? and is exactly the same as found from non-linear scattering 
in the scheme of the previous section. 

As the result is proportional to | £,|4 one must take all terms of the same 
order of magnitude into account and these come into the collision 
integral (right-hand side of equation (4.6)) from the product fst? 
x f°?) These give the Compton scattering and the interference of the 
Compton and the non-linear scattering. 

Although we derived the same results by another procedure, it is clear 
that there exist some questions about the singularities [w—(k-v)]~1 and 
€, 1, which indeed are not formal, but deep and physical. Still we find the 
stochastic heating and it is very clearly explained; it comes from sub- 
stitution Im [w—(k-v)]"! =—z 6(w—(K-v)) or a similar substitution 
for [o—@1—([k—k,]-v)]~, and this was done formally. Thus to clarify 
the physical meaning of the stochastic heating, one needs to treat the 
singularities correctly. 

It is known how this problem of the singularity [w—(k-v)]~? is treated 
in the linear theory. This is the Landau treatment that gives the Landau 
damping. It can be found as the limit of [o—(k-v)+ivr]-1 (where » is a 
small collision term) as » + 0. However, in the turbulent regime the 
energy in the fluctuations is much higher than at thermal equilibrium 
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and is, therefore, dominant. This means that the denominators must 
have a broadening dependent on the energy of the turbulence. Thus the 
whole Landau problem must be re-examined. This means that the nature 
of the Landau damping of plasmons is different from the usual linear 
damping. This also follows from the well-known fact that the linear 
Landau damping can be reversible and the plasmon Landau damping 
is irreversible. 

On the other hand, the presence of the energy of the turbulence in 
denominators that can be expected from this physical picture, means 
that to construct the corresponding theory it is necessary to sum the 
whole series in | £, |. 

This one can hope to do only because the turbulent collision frequency 
yr> is much smaller than the frequency of the plasmons, or at least 
than kv;. Thus we must use the small parameter 


elie ee (4.34) 


which is appropriate for weak turbulence. 


4.2. Turbulent Broadening of the Wave-particle Resonance 


To include the turbulent collisions it is necessary to change the per- 
turbation theory, and to start with an initial approximation that already 
is asum of the | E, |. This is possible to do as was shown by Rudakov and 
Tsytovich (1971), by a renormalisation of the equation for the stochas- 
tic distribution function, or by a renormalisation of the propagator 
{w—(k-v)}~). Let us write equation (4.7) for f° once more in the Fourier 
representation: 





he) stoch stoch k G) ~ 
—ife (kev) } fr +eE, (* Op 
= stoch Oe ) Al ] ge )y| 
--e| [ze (4. Fe oa Ey, ky Op / 
4 
x 4 dy, O(x—21— 22). (4.35) 


Let us extract then from the right-hand side of equation (4.35) the term 
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that is diagonal in f° and denote it by 


to —(ko0) + py fe = —e88( FL) sa mysc 
-e{ jes( i of) Coals oan am) 


Xx 6(%—%1—x2) d44 d4x2. (4.36) 
We introduce the operator g,(p) that is the inverse of the left-hand side 
operator of equation (4.36) 
8.(p) [o— (kev) + BP) fe" = fe, (4.37) 
and write equation (4.36) in the following form 


oe 1 Sra k fa sind 
yee = oe ade) (F an )+ +18. p)*(p) fer" 


+ gan f [egee(te 282) Caan. an) 


ki Op Am Nei ap 7) / 
x ieee dx; dx : (4.38) 


We insert this relation into the non-linear term on the right-hand side of 
equation (4.36) 


‘stoch 
— ilo — (kev) + ir,(p)] for" = ee of )+ ip) for 


k op 
k ky afres 
stoch y»stoch -Stoch ,-stoch\ 1 2 
iy ee ot ees 
vie | (BS meee otitay 2 yy (f. op) ()( Bp ) 


x 6(%—x1— 2) d'x, d4x_ 


: stoch{ ) 4 a oes 
~ie [ [a a “2 ) en sacorge : 


Op 
Stoc k 0 ‘Stoc 
~( Es 18 . ap) SP) >A DP) Se : ) O(x— 21 — x9) d4, d4x9 


3 k 0 stoc stoch ¢stoc 
rie | (E. sje wo (t +5) Les We h a toch p) 


SP ha Bes Be pyy (E; | Shaky Sosmagas pstoch py) ] 
xX 6(%—x1 —Hx2—x3) dtx, d4xy d4x3. (4.39) 
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In order that the operations that have been made have a meaning, it is 

necessary that the diagonal term, at least to first order, vanishes on the 

right-hand side of equation (4.39). This term comes from the last integral 

on the right-hand side of equation (4.39), that is, from 

{ Ee opps. dP (4. 40) 

We thus find an integral equation for the turbulent collision operator 
0p): . 
toc 

sup) fee = —iee SS [dn OY ise.) He. (4.4) 


Remember that g is related to » by equation (4.37). The new ee 
theory can now be constructed by using as a first approximation 


— fo —(k-v) +9,(p)) for" = — exse( aT). (4.42) 


From this we obtain 


stoch © pstoch k ee) 
fa a © Esc (EF), (4.43) 


This result differs from equation (4.13) in that [w—(k-v)]~+ is replaced 
by a propagator which involves the turbulent collisions. Equation (4.43) 
leads to a quasi-linear diffusion coefficient: 


D; = ie | an [ES Eee ae A &(p) (4.44) 


and a dielectric constant 
ed d3 (6) 
Sf ane be) (ee EY. (4.45) 


By this means we have included the effect of the turbulent collisions 
in the lowest order of perturbation theory. One can also work out the 
next approximation, that is, the equation corresponding to equation (4.22) 
(including equation (4.26) and (4.28)). The result is: 


a am a lero | 8, " |Ex, lg dx, 


+e le. val” [Eng [° [Eng ” 50-1 — x2) der 4x2. (4.46) 
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The quantity p® is derived from @® (see equation (4.28)) by replacing 
[w—(k-v)]"+ by g,(p). The quantity 2, ,, is derived from 2, ,, of equation 
(4.32) by replacing [w—(k -v)]~+ by 2,(P), except that the term correspond- 
ing to (4zi/k) oo)... _x,,x Must be omitted. [t may be shown that terms of 
this type are now included in the first approximation. 

We must now consider in more detail the equation for g. It is obvious 
that away from resonance the solution of equation (4.37) is the same as 
it was for the incorrect expansion procedure, that is, [«—(k-v)]~1. So we 
need only consider the solution near resonance. Introducing 7 = w—(k-v) 
and assuming that 

n <« Max {w(k), kur}, (4.47) 


one can consider g as a function of n, k, and p. The equation 
[wo — (kev) +i,(p)] 8p) = 1, (4.48) 


which follows from equation (4.37) is an operator equation and we find 
that the most important terms are the derivatives 6/07 that arise in 3, 
when the new variables are introduced. In this case g is diagonal and can 
be replaced by its eigenvalue, so that we obtain 


(a) 6) 
(n+. Ba D, "On ) x(n) = I, (4.49) 
where 
2 1 (ke ‘ 
D, a - fate & ¢ y | Ew ig ees k-k’ (p). (4.50) 


Because of condition (4.47) 7 can be neglected in equation (4.50) so that 


Dye io, (4.51) 
kek 
Pos us ue as [BP Bey, kk, Pky door. (4.52) 


It can be ascertained that the imaginary part of Do is small compared 
with the real part; we shall, therefore, assume Do to be real. 
If we now perform a Fourier transformation 


Lf ne 
gO, = oo Bek e”” dz, (4.53) 
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we obtain 


Og 0) 
Pee -+0Dog?, = —2ni 8(7). (4.54) 


The solution of equation (4.54) is 
gn k(p) = goexp(—3Dor'), 1 x 0, (4.55) 
while g,,, is determined by the condition 


(0), oo). 








Br, k 22k =—2ai. (4.56) 
There are two possibilities: 
(1) Do>0, 82,4 =0, go =—2ni; | nes 
(2) Dyed 2) 4 20, 95 = Dns. | 
In the first case 
gi e(P) = =i etree de, (4.58) 
In the second case ; 
Bk (P) = i{” elt +5 Dor? dr, (4.59) 


Introducing expression (4.58) into equation (4.52), one finds 


e 2 oo 
Do(k, p) = & { d4x, ey { dt exp & i{w,—(ky-v)} t 
1 0 


2 : Do(k—ki, py). (4.60) 


Since the integration in equation (4.60) is over a wide range of k and w 
one can approximate Do by 





ea fate (ke | Ex, [2 8(w —(K+0)). (4.61) 


This confirms that Do is real and positive. The second solution (4.57) 
leads to the same results, if one uses the substitution t— —t, as long as 
one takes the positive time direction to be always to correspond to in- 
creasing entropy. 


TSY 5 
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We can now also estimate the effective collision frequency. From 
equations (4.58) and (4.61) it is easily seen that 


Wa? k2y2 1/3 
1/3 pe' Te 
Pep ~ N~ Di? ~ (wae) 


nT kv 


kute 1/3 W 1/3 
“ on( 5] (si.) (4.62) 


This estimate is valid for electron—-plasmon interactions. A similar esti- 
mate can be made for ion—plasmon interactions. 

It must be emphasised that ».¢ cannot be expanded in the turbulent 
energy and that this comes about because of the importance of the turbu- 
lent collisions near resonance. Notice also that »,@ does decrease very 
rapidly with wave-number and that »,.@ is associated with the stochastic 
part of the distribution function. 

Let us now return to the problem of stochastic heating and the Landau 
damping of plasmons. The non-linear &, now describes such a damping, 
which in the zeroth approximation is identical with the linear Landau 
damping. If w #4 (k-+v), we can expand the equation for g, taking > as a 
perturbation. But g appears in ? and we now know how to treat it when 
resonance is possible. This is, in fact, possible for the beat frequencies 
and we can obtain an approximation by replacing g,_,,(p) in the integrals 
by 6(#—w,—([k—k,]-v)). We thus find the result obtained previously 
and can understand the physical meaning of the stochastic heating. 

We have not yet succeeded in avoiding all the singularities. Indeed, 
&_, In equation (4.46) has the same resonance as before, since it follows 
from equation (4.62) that the turbulent collisions produce only a small 
broadening and that the real part of &_, is almost the same as it was 
previously. 


4.3. Broadening of the Wave—Wave Interactions and Correlation 
Functions in a Turbulent Plasma 


We now suppose that the three-wave decay process is possible and 
that the singularities in equations (4.46) and (4.28) exist. The wave-particle 
collisions, even if they introduce an imaginary part into the non-linear 
é, do not remove the singularity. We must now renormalise the 
plasmon propagator (k2—we°)—1, 
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This problem was considered by Makhanov and Tsytovich (1970). 
Here we give only the main ideas. It is obvious that one cannot use 
equation (4.27) to express the average of the three fields by an average 
of four fields. 


Two steps are essential. The first is the introduction of a non-linear 
(renormalised) dielectric constant e? = ~ | dtd, |E,,|?. The second is 


to split the stochastic field into two parts: the high-frequency part E™, 
with frequencies of the order of the mean plasma frequency, and the 
low-frequency part E’ with frequencies of the order of the frequency 
of the plasmon-plasmon interactions (that is, the frequency corresponding 
to the correlation time of the fluctuations). We rewrite equation (4.20) 
as follows: 


4n ff ~¢ : 
(é, 4 Aaa) Oneal = : A ( Exochestoch = Ene) 


eS i} 
° UH, Xa) 2 m9 
ik 


4 4 4n f ~0 toch ;,stoch ,.stoch 
x 6(% — 1 — 2) d 41 d Met | Deve Rae one 1 ae 1 Da 


stoch / ;stoch ,,stoch stoch ,.stoch ;~stoch\, 
—£,, (Eur Ex ) — (En Ex, E, / 


xs 


x 8(« —x1 —x2—H3) d’x1 d’xe d’x3— i 2a Pes Sree Dae Soa) 
x 6(%—x1 — Hg — 3) d'x, d'x9 a3. (4.63) 


The ¢@ include the turbulent collisions considered in the previous section 
(although this is not crucial for the considerations given below). Thus by 
averaging equation (4.63) multiplied by E, we get a result similar to 
that of equation (4.22) with &,+é2" on the left-hand side instead of 
&, and the term e*' |E, |? on the right-hand side. By definition, the last 
term must vanish to take into account all the terms diagonal in |E,|? 


on the left-hand side. The average of the three turbulent fields can be 
written in the form 


(epee) = (etete) + (eteles) + (etetel): 


2 
Eso. EM 4+ EY, (4.64) 


The E,,-field in equation (4.22) is always a high-frequency field. Thus 
one needs to express the low-frequency field in terms of the high-frequency 
5¢* 
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field E4. This can be done if one takes into account that a low frequency 
can come from the difference of two high frequencies and thus 


EL = 1 sis ae { eee (EES — (CERES ) 6(x— 2, — x2) d‘x, dsp. 
dy 


(4.65) 


We can also of course write down the next terms of equation (4.62), 
but this will not help in the calculation of expression (4.64). Thus, by 
using equation (4.65) in the last two terms of equation (4.64) and also in 
equation (4.25) one finds from them one of the terms diagonal in |E,|?; 
the other one comes from 6®), The requirement that these diagonal 
terms do not occur on the right-hand side of the equation which we 
derived leads to the relation 

v _ 4ai ™(3) _Am(4n)’ Dee Hoy Or, x, =x, (4.66) 
k|k— —ky| By wy t Boy 





This equation together with the definition of &-" gives a non-linear 
integral equation for e™', Because the denominator Se ae occurs 
in an integral in the induced decay process, one can sometimes neglect 
é™!. and put approximately 


Im — Eres ~ — 7 O(Ex). (4.67) 


&, +e 


To calculate (EFESES) one needs (i) to write down the equation 
. 63) for E¥, taking on the cae hand side only the terms linear in 
L (ii) multiply the result by EE," H; and (iii) average over a statistical 
oe Then, by using equation (4.65) and approximating the five-field 
average terms by a product of the averages of three and two fields, one 
finds a non-linear equation: 


Se bot 8ai Es ‘ 
(Bet Fi") (Ev Ee Eas) = pr Be) ~nyy—ny [Eng ? [Eng 5G! + ea + 222) 
+m Help 
. |B, |’ so HQ, My , (ESERE, 2 6(x’ +1 — x9 —%3) 


x d'x5 d’ug+1 2 ai, (4.68) 
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where 
~(2) ~(2) 
~(3) 820s, + x5, Hos x0x, My Xa tHy (4 69) 


~X, H1) Ht2,%3 = 0x Hy, %2, XB ~n. |. . 
, |kot+ks| (é Eeatug t Enetxg) 


x 


One can solve this equation by a perturbation method using only the 
first term of the right-hand side of equation (4.68) in first approximation 
and one finds that the other terms give small corrections of the order of 
W/nT < |. Thus the equation for the correlation function takes the 
form 


~(2) 12 
(42) |E,2 = 2 or [Be et 


X 8(%— 414 —H2) d4xy dx. (4.70) 


This has (i) a soluble form as |E, |? => 0 and there are no secular terms, 
(ii) gives the same expression for the balance equation, and (iii) gives 
also the same tail of the correlation function as without €""- 

Thus we have seen that the turbulent collisions are essential in the 
treatment of the properties of a turbulent plasma and the appropriate 
treatment of them can be found by summing an infinite series of some 
of the terms which are essential near resonances. The turbulent broaden- 
ing can be neglected only for such average values as the diffusion coeffi- 


cients or the expression | |E, |? dw. In the next chapter, we shall show 


through the example of ion-sound turbulence how essential this turbulent 
broadening is, which is due both to wave-particle interactions and to 
wave-—wave interactions. 


5. The Spectrum and Correlation Functions 
of lon-sound Turbulence 


5.1. Introduction 


Ton-sound turbulence can be excited by different processes. For 
example, it is well known that such turbulence can be excited in shock 
waves propagating perpendicular to an external magnetic field or arising 
from an external electric field along the magnetic field lines. In both cases, 
excitation arises if the mean drift velocity uw of the electrons is larger 
than the ion-sound velocity », and T, >> T;. A detailed numerical calcula- 
tion of the growth-rate of this instability was given by Stringer (1964). 
The essential behaviour of this kind of instability is that, if u > v,, the 
growth-rate exceeds the Landau damping for each frequency. Dissipation 
possibly exists only for the waves that propagate outside the Cherenkov 
cone or, in the case when u >> v,, in the direction opposite to the applied 
electric field. This means that an exactly stationary turbulence can be 
created, if there is energy transfer in the angles of the turbulent waves. 
Another example is the excitation by non-linear creation of waves that 
can be due to the decay of waves with higher frequencies (for instance, 
Langmuir waves). [In this case, the excitation can exist in a narrow 
interval of frequencies so that the lowering of the frequencies of the ion- 
sound wave can lead the system to a stationary turbulent state. 

An essential question is what kind of non-linear interactions can be 
involved in such an energy transfer. Usually, ion-sound turbulence is 
accompanied by an anomalous resistivity of the plasma, due to the scatter- 
ing of the electrons by the turbulent ion-sound waves. This is usually 
observed and measured experimentally. The value of the anomalous 
resistivity depends, of course, on the nature of the non-linear energy 
transfer. 


In principle, one can suppose that there can exist the three following 
non-linear processes: 


e+s = e’+s', (5.1) 
i+s + i’+s’, (5.2) 


st sts". (5.3) 
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The first two processes are induced scattering processes and the last 
one is a decay process. The process (5.2) was first considered by Kadom- 
tsev (1964) in a partially ionised plasma. The condition that the ionisation 
of the plasma is very small seems to be a very essential point in the 
Kadomtsev theory, that is, the binary collisions must play a significant 
role. Indeed, as a result of the presence of a high number of collisions, a 
region where damping exists is created at not very low frequencies. 
Energy transfer to this region due to a lowering of the frequencies is 
possible. Still, in the Kadomtsev treatment as given in his book (1964) 
it is supposed that the waves mostly go into a direction near to that of 
the applied field because, for example, the plasma is confined in a long 
tube, and the waves that move at larger angles strike the wall of the tube 
and are, therefore, damped. The turbulent energy is not very high 
because of the binary collisions. These two restrictions (the collisions and 
a long tube) are very essential. Nevertheless, some authors tried later 
to apply the same equations to a fully ionised plasma, or to the case 
where the change in the angle of the wave can be rather large. Indeed, 
from the probabilities derived earlier we find the following non-linear 
growth-rate “due to the scattering of ions”: 








a, ____Im etl. 
rn 
OO |u=uy 
yn 8 Ti Of Ce kn)(UK A bal Th A ki) 3, 
=k ak T, fas Series ~kkinT. mx oa Wy, 20k Us, 
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x, = ED, = Vein) (5.4) 
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One can see from equation (5.4) that if there exist waves propagating 
in a direction close to that of the applied electric field, waves with fre- 
quencies sufficiently lower than the original ones are growing in any 
direction including the direction opposite to that of the field. Thus, the 
distribution of turbulent waves must be more or less isotropic. This seems 
to contradict the known experimental results (Janéarik and Hamberger, 
1970; Daughney et al/., 1970). On the other hand, the interaction (5.4) 
seems to be rather weak. That means that to compensate the high 
Landau damping one needs a high level of turbulent energy. Indeed, if 
one takes into account that the linear decrement is of the order of 
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V($2m,/m;) and requires that it balances expression (5.4), one finds the 
following rough estimate for the turbulent energy: 


W any Me T 2| max 

7m 4 (7) ain Ea (5.5) 
which shows that W seems to be too high (as we assumed that W/nT, « 1). 
Estimating W from balancing the non-linear interaction with the linear 
growth-rate of the current driven turbulence—which is of the order of 
V(s7)-(m,/m;) (u/v,) o—and putting it in the quasi-linear equation makes 
it possible to estimate the turbulent collision frequency of the electrons 


and, therefore, the turbulent conductivity. This was done by Sagdeev 
(1967), who found that 
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th, (5.6) 


where «, when calculated from equation (5.4), is of the order of unity 
or even less. To obtain agreement with the observed thickness of shocks 
it is necessary to put in equation (5.6) an experimental value for « of 
the order of 102. This means that the non-linear interaction must be much 
more effective, if one believes that the observed anomalous resistivity is 
due to the development of the ion-sound instability. 

The problem became much more sophisticated after Drummond’s 
suggestion that the most important non-linear interactions may be due 
to the scattering of ion-sound waves by electrons, that is, to the process 
(5.1). Indeed the rough estimates by Sisonenko and Stepanov (1969) and 
by Krall and Book (1969) show that such an interaction produces a 
balance for the linear growth-rate, even if the level of the turbulent energy 
is very small, 

W me u 

WE. ae (5.7) 
This also contradicts the observations that show the turbulent energy 
level to be high (Janéarik and Hamberger, 1970; Daughney et al., 1970). 

The interaction (5.3) is usually neglected because the ion-sound waves 
do not satisfy the conservation laws 


w(k) = w(ki1)+ (ke), (5.8) 
k=ky+kg. (5.9) 
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If this process is allowed, it dominates the processes (5.1) and (5.2), 
because the scattering arises as a tail of the decay process. 

To finish this introduction, it is necessary to mention the work of 
Rudakov and Korablev (1966) and of Kovrizhnyk (1967) who considered 
the quasi-linear development of ion-sound turbulence. The argument 
that the quasi-linear interaction may be a most important one can be 
seen as follows. The non-linear interaction (5.4) is rather weak. Therefore, 
one may assume that the quasi-linear change in the spectrum is unimpor- 
tant. On the other hand, the quasi-linear development can also change 
the angular distribution or, in other words, bring the oscillations from 
the excitation region to the damping region. However, the best known 
quasi-linear effect is to form a kind of plateau. In this case, it means 
that the growth-rate comes near to the threshold, that is, », ~ 0. This 
means that wis of the order of v,. On the other hand, as shown by the 
above authors, the spectrum of the turbulent frequencies must be very 
narrow. The two statements (u ~ v, and narrow spectrum) seem to 
contradict the above-mentioned observations (broad spectrum and 
essential difference between u and 2,). 


5.2. The Influence of the Turbulent Collisions 


Now let us consider from the general point of view of the turbulent 
collisions the interactions with ion-sound waves, and illustrate the essen- 
tial physical change of the picture of the turbulent state, if these are taken 
into account. First of all, the interaction of ion-sound with the electrons 
is taken into account in the quasi-linear interaction. We cannot expect 
that the next term in the expansion, that is, the non-linear interaction, 
can be expanded in terms of the turbulent energy, because the turbulent 
collisions were already taken into account, and the result cannot be 
expanded in terms of W. In other words, only the first rough approxima- 
tion corresponds to the quasi-linear one as was shown above, and the 
exact expression cannot be expanded and, therefore, the corrections to 
this rough approximation can also not be expanded. This means that 
the broadening of the resonance must be supposed to depress the non- 
linear interaction, that is, the difference between the rough and the exact 
solution as well as higher-order effects. Indeed, in any consistent theory, 
in which there is a small parameter, the next order (perturbation effect) 
must be small, so that one can expect the effect of process (5.1) to be 
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negligible. This was indeed shown by Rudakov and Tsytovich (1971) to 
be the case. 

The next question is the possibility of the process (5.3). It is forbidden, 
but not strictly so. Indeed, if one ncglects the curvature of the ton-sound 
wave branch and supposes that 


ws, = kvy, (5.10) 


then for three waves propagating in the same direction the fulfilment of 
the conservation law (5.9) means the fulfilment of the conservation law 
(5.8). However, if there is a curvature, this is no longer possible. What 
is the effect of the turbulent collisions associated with the wave—wave 
interactions? If the spreading-out of the dispersion curve due to these 
turbulent collisions is larger than the effect of the curvature, one must 
say that such a process is allowed. Because the probability of this process 
is large, one can suppose that the non-linear interactions become much 
stronger. This is the effect needed to obtain agreement between the theory 
and the observed thickness of the shock waves. 


5.3. The Electron Ion-sound Non-linear Interactions 


It is easy to see that the large cross-section for electron ion-sound 
interactions found by Krall and Book (1969) and by Sisonenko and 
Stepanov (1969) can exist for a very small level of the turbulent energy, 
much smaller than that for which this non-linear interaction can balance 
the linear growth-rate. Indeed, the perturbation theory expansion used 
in these papers can only be supposed to have any meaning if the particle 
oscillation velocity in the turbulent field ». is much smaller than the 
mean particle velocity. In the case of ion-sound, one can write 


ef stoch eFstoch 


~~ aloo eo) ~ “ny” a 


because the Cherenkov condition means that (k-v) > @,, while the lowest 
(K-v) is of the order of w,. Therefore, 
i 2 EAE MOEN ae VE _ in W 
UTe . mkwjv?,. = wo; nT. ~ me. nT, , 
Here w, is of the order of w,;. We find thus that even from these argu- 
ments one can see that the limit (5.7) which is necessary to balance the 





(5.12) 
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linear growth-rate cannot be reached. However, one can see that the non- 
linear interactions are lowered by the turbulent collisions even for much 
smaller W. 

Indeed, the theory of turbulent collisions developed in the previous 
section is a theory that uses the small parameter 

ae e%( Estoch)2 z w2.W 
me) oF. ~ ATO? 

If one uses the estimate for k ~ c,¢/v7_ (the most important k in the ion- 
sound spectrum) one finds that 


W \18 
=(.7,) (5.14) 


Therefore, we have for weak turbulence that 2 « 1. As we shall see now, 
the corrections to the first approximations, that is, to the quasi-linear 
result, are larger than the result (5.14). Thus, the second approxima- 
tion need not be considered. Therefore, the non-linear interactions of 
ion-sound waves with electrons come from the difference between the 
exact expression of the first approximation, taking into account the 
turbulent collisions and the approximate one when Im g is replaced by 
—n6(w,—(k-0)). We thus find that 


yas _ P(E) bm 8x80 0 be 
aa Ile > Hoe (h0) rr | 


where 7, is an approximate growth-rate. Using the expression for g one 
finds, for example, for a Maxwellian distribution 


Ch rs 2 
Pee oc a exp (+ a a 
0 


(5.13) 


(5.15) 











Vk wu 
WT , ot at? | wt at 
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where _ Dov 
= (5.17) 


In t there exist two essential cut-offs, namely t = 1, ~ ykv7,/m and 
T= te © y"3 fol3, If tp >> t1, one can put the exponential equal to unity 
and one gets 





ee aL (5.18) 
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This is approximately the result of the rough estimates of Krall and Book 
(1969) and of Sisonenko and Stepanov (1969). 


If t2 < 7, one finds 
n.d. 1/2 
Vere x (wi) (5.19) 


These corrections, indeed, are larger than expression (5.14). The result 
(5.19) shows that the non-linear corrections due to the electrons never 
can reach the linear effect, if the turbulence is weak, W/nT, « 1. One 
can now estimate the value of W for which t; = te, and one finds that 


it corresponds to 
2 
ee ee a (5.20) 


This is much smaller than the estimate found from equation (5.12): 


pt ates (5.21) 








Fic. 16. 


The ratio of the non-linear correction y®"" to the quasi-linear growth-rate 

y,, as function of the ratio of the turbulent energy W to the thermal par- 

ticle energy nT,. The solid curve takes into account the turbulent collision 
broadening. 
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16). If one linearly extrapolates the results found for very low W, that is, 
Ww m.\? 
a< (Fr) (5.22) 
to the value when this result is not valid, one gets the wrong answer 
that the scattering by electrons is effective, that is, v2" /y, ~ 1, when 
W ~~ me 
nT. om; 
This illustrates the significant role of the turbulent collisions. 


(5.23) 





5.4. The Wave—Wave Interaction of Ion-sound Oscillations 


This kind of interaction was supposed to play a significant role in the 
formation of ion-sound turbulence by the author (Tsytovich, 1971) in 
a paper written during his stay in Oxford and Culham; this was the result 
of discussions on the experimental work done in the Culham Laboratory 
by Janéarik and Hamberger (1970) and by Daughney ef al. (1970). We 
describe here only the qualitative picture of such an interaction. One can 
see that, if the level of the turbulent energy is high enough, the process 
(5.3) is allowed, if the correlation width is of the order 


> (Se). (5.24) 
@pi 
This means that the necessary conditions can be satisfied if one goes to 
lower frequencies. If such an interaction is involved, the waves created 
do not change their angle during the non-linear energy transfer. The 
interaction (5.3) conserves the plasma energy and does not change the 
angle. The nature of this interaction is to spread the energy over a larger 
q-interval without changing the total energy. Therefore, it creates a broad 
spectrum. But the condition (5.24) means that in frequency space there 
exists a wall, because for some frequency greater than the maximum 
frequency defined by equation (5.24) the process (5.3) is forbidden. Thus 
the broadening can only be the process of the lowering of the frequency. 
One can see that the spectrum created in such a process is in some sense 
similar to the ion-sound ion interactions considered by Kadomtsev: 


1 8role 
W, eer  peP(A), (5.25) 
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where the quasi-linear growth-rate is written as wye, 


_ (k-E) 
oon kE ° 





(5.26) 


so that y, describes the non-dimensional growth-rate and P(A) is a loga- 
rithmic-type function varying slowly with A = ,,,/@. The difference 
from the Kadomtsev interaction, apart from this logarithmic function, 
is a greater efficiency of the interaction (5.3) and, therefore, a lower value 
of W needed to balance the growth-rate. However, the correlation width 
as an effective turbulent collision described in the previous section in- 
creases quickly as the frequency in the spectrum decreases. This gives an 
increase in Mw and, therefore, the possibility of interactions of waves 
with different directions of propagation. Indeed, this is allowed by the 
conservation laws for angles: 


Ap =~", (5.27) 


One can then estimate the time needed to transfer the turbulent energy 
by a multiple change in the angle, each change having the value (5.27), 
to the final result of an angle change of the order of unity. 

This time can be compared with the time needed to propagate the wave 
energy into lower frequencies and one finds that the highest frequency 
does not much differ from @,; and the lowest frequency into which the 
energy can go before the angle changes become of the order unity, is 
of the order of 3 ¥/(m,/im;)@,;. When the angle is changed in the current- 
driven instability, it means that the energy goes to the dissipation region 
and can be absorbed. Of course, in this process the ions are involved and 
the ions or electrons are heated, but mostly the electrons.The tail of the 
decay interaction is the ion-sound scattering by ions. It is possible that a 
multiple change in angle due to correlational broadening has a greater 
probability than the change due to one step of ion-sound scattering by ions. 


5.5. The Anomalous Resistivity of the Plasma 


Even if the whole spectrum is not stationary, the anomalous resistivity 
can obtain a stationary value, because it is determined by the highest 
frequency in the spectrum. Thus, to have a stationary, or more precisely 
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quasi-stationary, resistivity it is necessary only to form a spectrum near 
the maximum frequency. Then one can write down the quasi-linear equa- 
tion that describes in first approximation only an elastic scattering of 
particles by waves or only the change of the angular dependence of the 
distribution function of the electrons. This follows from @/kv7.<« 1, which 
is found for ion-sound waves. 

If x is (p-E)/pE one finds that the quasi-linear equation is of the form 





V(l-x) » Off 
0% _ 2m Ve O dt | dw Gere (5.28) 
* op a Meds Ox —~V(~x?) Vf (1 — x? &) . 
where 
W= 2a { W.,2 dE do. (5.29) 


We have split the electron distribution function f* into an isotropic 
part fo°® and an x-dependent (anisotropic) part f{*%: 


free = freq free, fre free, (5.30) 


On the other hand, W,,,. is proportional to yz, which depends also on 
the function f{*: 


On va-#) dx, 
a ae ey as [MOM Gap 63D 





The self-consistent solution can be written in the form 


Sie = g(x) =. ee Uss (5.32) 
where 
29) VF 99, s=y(l-x), «= —___., 6.33) 


émarsonu | Wadd 
while @(s) satisfies the equation 


sf} e(n) dn ea (s+n)? nt 





cGy (5.34) 


The numerical solution of this equation is shown in Fig. 17. 


os) J, Vn) Gon 


72 


THE THEORY OF PLASMA TURBULENCE 


0:8 
0-6 
gs 
a 04 
02 
: s 
0 0:2 0-4 06 08 1:0 
Fic. 17. 


The result of the numerical solution of the non-linear equation for 9(s). 


From o(s) one can easily find the angular distribution of the turbulent 
spectrum to be proportional to A(é), where 


aE) = £ i) : eee (5.35) 


The result of a numerical computation of A(€) is shown in Fig. 18. To 
find the conductivity it is necessary to find f® from equation (5.32). 
One sees that then 


(5.36) 


(5.37) 
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The angular dependence of the ion-sound spectrum, 


equation (5.36) can be written in the form 


ner 
c= ae (5.38) 
where 
1002, 


This value seems to be in agreement with the known data on the 
thickness of shocks and the measurement of the anomalous conductivity 
in an electric field. We must bearin mind that w,,, is determined by the 
correlation time or by Aw (see equation (5.24) ) and that thus the 
conductivity depends on Aw as follows from equation (5.39).t 


t The spectrum calculated here extends only to about w,,/2. To find the con- 
ductivity, the spectrum mus tbe extended to w,, (after that it drops off rapidly). 
There is no experimental evidence that the spectrum beyond w,),/2 produces the 
major part of the conductivity, and it is therefore unlikely that the spectrum will 
have a large peak between ,):/2 and @,). 


TSY 6 


6. The Spectrum and Correlation Functions 
of Langmuir Turbulence 


Tue Langmuir plasma oscillations are the most typical collective plasma 
motion. The frequency of these oscillations in a sufficiently weak magnetic 
field, when @y, < Wye, is approximately constant, with a weak wavelength 
variation which is due to the thermal particle motion: 


dp toes. 2, (6.1) 


There exist different types of mechanisms to excite such a motion. One 
of the best known is excitation by electron beams. This excitation is 
due to a population inversion of the particles if the beam is present, and 
the possibility of Cherenkov resonance with the Langmuir waves. The 
beam does not greatly change the dispersion law (6.1) for the waves 
provided the density of the beam, 7, is sufficiently small, and the spread 
in the particle velocity Av sufficiently high: 


my (F)- (6.2) 


Thus, the presence of such a beam changes only the imaginary part of 
@,, that is, such waves are excited. It is not always necessary to have a well 
directed beam. Any anisotropy of fast particles with v >> v7, will produce 
instability of the Langmuir waves. Analogously to the case of ion-sound 
waves it is possible to have also an excitation of Langmuir waves by the 
non-linear decay of waves with higher frequencies. Such waves are the 
normal transverse waves. Thus one can say that a high frequency field 
and laser beams can excite these turbulent motions. 

Indeed, these two mechanisms of excitation of Langmuir waves by 
the beam and by transverse waves are similar, as one can see from the 
diagrams of Fig. 19, which describe the Cherenkov and decay resonance. 
The mathematical description of these two processes of turbulent 
excitation is also similar. The essential difference between ion-sound 
jnstability and Langmuir instability is that the growth-rate of the latter 
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Emission of Langmuir waves by particles (a) and by transverse waves (b). 


is usually restricted to a particular k-interval, for example, near the phase 
velocity of the beam. It is obvious that these velocities are less than the 
velocity of light. In the case of the non-linear excitation of Langmuir 
waves by a “beam” of transverse waves one finds that the excitation 
is forbidden for waves with phase velocities larger than the velocity of 
light. This results from the fact that the group velocity of the transverse 
waves is less than the light velocity. The higher the beam-particle velocity, 
the more important the non-linear interactions that can stabilise or lower 
the growth-rate of the beam-—plasma interaction. The excitation of 
Langmuir waves with relatively small phase velocities is thus the most 
efficient process. For the non-linear excitation the most efficient process 
is the excitation of Langmuir waves with a wavelength of the order of 
magnitude of the wavelength of the transverse wave which produces 
the excitation, so that, when w > w,,(c/v7,), the excitation by non-linear 
processes becomes inefficient. The stabilisation of this kind of excitation 
by transverse waves due to the non-linear Langmuir wave interactions 
is more efficient than for the beam—plasma interaction, because the 
growth-rate of the non-linear excitation is not very large. 

There are the following three possible non-linear processes that can 
create the turbulence spectrum: 


etl=ze'+l, (6.3) 
itl =i'+I, (6.4) 
lah =a (6.5) 


Here | denotes a Langmuir wave and e and i the electrons and ions of 
the cold plasma, respectively. One can include also the process of the 
scattering of Langmuir waves by the beam electrons in the case of a 
beam-plasma instability. However, this comes as a non-linear correction 
to the quasi-linear growth-rate, which itself is already small. The reason 
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why the process (6.3), that is, the scattering by cold electrons, must be 
considered is that the cold electrons are not in resonance with the Lang- 
muir waves which are excited. 


6.1. The Spectrum of Small-scale Langmuir Turbulence 


Estimates show that for high k-numbers (or small dimensions) the 
effectiveness of the induced scattering by electrons decreases, when k 
decreases, as k*, and is most effective for the largest possible k. This 
value of k cannot be larger than 1 /Ap, where Ap is the Debye length. The 
scattering by ions decreases, when & increases, as k~2, and, therefore, 
one can estimate where the electron scattering is the dominant process 
and finds 


] ] Me V5 
aS > k > ke = an (sn . (6.6) 


The plasmon-plasmon collisions (6.5) are negligible in the range (6.6). 
Only in a plasma with sufficiently heavy ions is this interval sufficiently 
broad to be considered as an important region. Also it is necessary that 
the excitation of turbulence occurs very near k ~ 1/Ap in order to have 
a region where one can neglect the excitation and to be justified to assume 
that the spectrum is due only to the balance by non-linear interactions 
with the electrons. Nevertheless, one can work out the possible turbulent 
spectrum in the range (6.6) in order to have a complete picture of the 
Langmuir turbulent spectrum in a plasma. If one supposed that, for 
example, the distribution of Langmuir waves is isotropic, it is easy to 
find from the probability for the process (6.3) the non-linear balance 
equation 


OW, | oo eae Gee ; 
ay = airs | [ Wi, dk 9 (B-R) (- + K) 
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B= Sure V (2m) | (6.7) 
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After some transformations and differentiations equation (6.7) becomes 


BW, 21 @W, 46 dW,  10W, 


Ge 6 Ge le dee ae se 
There are three solutions of equation (6.8): 
_ const | a, _ stvVI7 
Wie Dis US (6.9) 


Only one of them gives non-divergent integrals in equation (6.7) and, 
therefore, it is the only solution of equation (6.8): 


p= 4, (6.10) 


In many cases, the beam growth-rate is not zero in the interval consid- 
ered. In those cases, it is necessary to find the quasi-linear change in the 
growth-rate which is due to the quasi-linear beam relaxation. This 
means that in the case of an infinite plasma and an infinite beam inter- 
action the power input is not constant and one must take into account 
also the change of the turbulent energy that is due to the change in input, 


6.2. The Spectrum of Langmuir Turbulence in the Intermediate-scale 
Region 


The intermediate-scale region is defined by the inequalities 


k,<xk<k,,, (6.11) 


tak I Me 
ke =yory Vi (6.12) 


In this region the scattering by ions is much more efficient than the 
scattering by electrons. On the other hand, the process (6.5) need 
not usually be included. The exact non-linear interaction with ions is 
described by the equation 


ONK N, 3 Opel eT; 


at * 8 ngmeor (1+ (7.7? 
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(6.13) 


where 
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which can be found from the probability for ion-plasma scattering. One 
can see from this equation that the non-linear growth-rate increases 
with the frequency difference of the interacting waves until this difference 
is lessthan |k —ky| v7;. If |k —k1| is of the order of |k| (which can occur 
even when |k,|isequal to || in the case when the angle between & and ky 
is not small) one finds the maximum k,—k = Ak for which the transfer 


is most probable 
Ak _!t me\ kvte 
ko a VR) “Opes (ly) 


Thus k/k <1 when k<k,. If then one considers 4k = k, to be a 
physically infinitely small value in the interval k >> k,, one can find a 
differential equation for the energy transfer. 

Indeed, supposing that the turbulence is isotropic and using the 
approximate description of the function 


OF ep ae ae (6.15) 
~ Tea eg, °° |~ aaa] = 

one finds an equation 

OW, _ aw, 

or ae VWrtaoW, BE? (6.16) 
where 

~ TO, ete. 
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Such an equation can be applied both in the region where the excitation 
exists and in the absence of excitation and damping, when y, = 0. If 
¥, = 0 and approximately constant in the small k-interval 4k one finds 
the linear growth of W, as shown in Fig. 20. 

Outside Ak equation (6.16) gives W, = const. The value of this constant 
is determined by the value reached in the interval Ak. Calculating the 
power of the turbulent generation 


Q= | mare dk, (6.18) 


where the integral in equation (6.18) is only over the k-interval where 
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The spectrum of Langmuir turbulence created by excitation in the range 
Ak and scattering by ions. 


y, # 0, that is, over Jk, one can express the constant in the spectrum of 


Fig. 20 in terms of Q: 
_ 4/(22 
W, = (z= ) (6.19) 


6.3. The Spectrum of Langmuir Turbulence in the Large-scale Region 


The large-scale interval is determined by the inequality 
k<k,. (6.20) 
In this region the turbulent energy can be transferred directly to very 
small k-values. As all the plasmons are pushed into a very small phase- 
volume, the plasmon-plasmon scattering (6.5) becomes very important. 
This interaction is accompanied by the scattering of plasmons by ions. 
Therefore, in the formation of the turbulent spectrum both processes (6.4) 
and (6.5) are involved. It must be mentioned that the probability for 
the process (6.5) increases steeply when k < ks. In the region k < kya 
power-type solution for the turbulence spectrum was found by Pikel’ner 
and Tsytovich (1968): 
Wi = oo, (6.21) 
where » satisfies the equation: 
] 1 1 1 
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where F is the hypergeometric function. If Q is sufficiently large, one can 
neglect /°; (which describes the contribution from the process (6.4)), 
and the solution of equation (6.22) ts: 


yp x 2-84, (6.23) 


This solution is found as the intersection when the right-hand side and 
the left-hand side of equation (6.22) are plotted as functions of » as 
shown in Fig. 2!. From Fig. 2! one can see that an appreciable change 
in /* does not greatly change ». If J; is much larger than Io, one sees 
from Fig. 21 that » will tend to the value 4. 
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The left-hand side (curve !) and right-hand side (curve 2) of equation 
2! 


The analytical result was checked by computer analysis of the same 
equation. In the asymptotic region good agreement with the analytical 
solution is found and this is illustrated in Fig. 21. 
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(a) The spectrum of Langmuir turbulence in the large-scale region. The 

solid curve represents the results of numerical calculations, and the 

dashed curve the approximate solution put on the computer. (b) Com- 

parison of the analytical asymptotic solution (dashed curve) and the 

numerical solution including the energy-containing region. The coinci- 

dence between the analytical and the numerical solution starts at 
k~(4to 5)ky; AQYU™ is the numerical valuc of ky. 


For the computed value of the turbulent energy the analytical solution 
gives the experimental index 
vp = 3-86, (6.24) 


while the result of the computer calculations gives 


y = 3:89. (6.25) 


The calculations give the whole curve of the turbulent spectrum, including 
the energy-containing region k ~ ko. The value ko must be much smaller 
than k,. The analytical solution makes it possible to estimate ko by 
joining the spectrum in the different k-intervals. In this way one finds 
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that 





Q 


This result was used for the approximating initial curve in the computer 
calculation and the resulting curve of the turbulent spectrum has a 
maximum very close to that of the initial curve (see Fig. 22). 


1\ 1720-1) 212 
pee k,/ Q ) sone a (6.26) 


6.4. The Spectrum of Langmuir Turbulence in a Non-isothermal Plasma 


If the temperature of electrons is much higher than the ion temperature, 


T. > T;, (6.27) 
a new non-linear process is involved in the formation of the spectrum: 
l= I'+s. (6.28) 
This process is allowed only when 
k>k, (6.29) 
where 
ki = ke Va (6.30) 


The decay process (6.28) dominates over both the plasmon-ion scattering 
and the electron-plasmon scattering, that is, both the processes (6.3) and 
(6.4). This is a result of the fact that the scattering processes always are 
the tails of the resonance decay process, if such a process is allowed. 

Therefore, the region W, ~ k~5/? disappears in the case when T, >> T;. 
If k >> ks, the energy transfer in the process (6.28) again is described by 
a differential equation, while 


—— : (6.31) 
The equation that describes the energy transfer due to the decay (6.28) 


is the same as equation (6.16) with a change of the constant « toa’, given 
by the relation 


r_ T. \? 
4 =o(1+7") F (6.32) 
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This means that for the same power Q the constant value on the turbulent 
spectrum is smaller by a factor 2 in a non-isothermal plasma with 
T, >> T, than in an isothermal plasma, where T, = T;. 
In the region 
k<k<k (6.33) 


the decay is forbidden and the scattering by ions is the most important. 
The turbulent energy in this interval thus increases steeply, because of the 
difference of a factor (T,/T,)? between the constants « and a’. 

In the region where k < k, the same solution as that found for an 
isothermal plasma occurs. The value of ko which describes the position 
of the maximum of the spectrum is in the case when T, > T7; larger than 
for the case where T, = 7;. This is illustrated by Fig. 23. 
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Sketch of Langmuir turbulence spectra. Curves 1 and 2 are given for the 
same Q, but curve 1 corresponds to T, = 7; and curve 2 to T, > T;. 


6.5. The Radiative Type of Langmuir Turbulence Spectra 


In the previous considerations the turbulent energy in the energy- 
containing region was damped by ordinary damping due to binary 
collisions with the plasma waves. This can be the case when the plasma 
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is (i) optically thin for the radiation processes and they do not play any 
essential role in energy transfer or (ii) optically thick where the radiation 
emitted by turbulent motion is reabsorbed and converted back into 
plasma waves, and, therefore, the energy in the transverse waves is 
negligible and they do not change the turbulent spectrum. 

The optical thickness L, of the plasma for the non-linear emission 


I+) =t (6.34) 
can be estimated as follows: 
=. “/ Me OpcltoT 
Le = yom / ln ‘ ) (6.35) 


If the plasma characteristic dimension L is much larger than L,, the 
process (6.34) is unimportant in the formation of Langmuir turbulence. 
If L< Ls and 0 > QO, where 


cc 


; Me. 
iad Reo Se Wpellole 6 
mM, UTe 


1 


2 

€ 
aa ; (6.36) 
the radiation losses dominate over the collision damping. Thus, under 
such conditions most of the energy that is put into turbulent motion goes 
into the radiation, and the turbulence is of a radiative type. The change 
in the spectrum is different in the case where the radiation losses are 
able to damp the turbulent motion before they reach the maximum of 
the spectrum in Fig. 22 from the case where this damping becomes 
essential when the turbulent oscillations are able to reach this maximum. 
In the first case, the spectrum is flat: W, = const for k > k,, and 
exactly the same as if the radiation losses were absent (that is, the expres- 
sion (6.19)) up to the maximum kya, for which the radiation losses 
give an exponential cut-off of the spectrum (see Fig. 24; Tsytovich, 1969a) 
The ratio ka, to k, is of the order of 


Kei aS 2 m, Ute T : 
L & 10 ee (1+ 7°) : (6.37) 


The situation shown in Fig. 24 occurs when Kyay/ky >> 1. 1f Kmax/Ke <1, 
the energy reaches a maximum and the exponential cut-off disappears. 
A maximum is formed which depends weakly on Q. However, the energy 
loss is due to the emission. One finds that the damping due to emission 
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The Langmuir turbulence spectrum for the case when the turbulent energy 
is dissipated into radiation (T, = 7)). 


is described by the equation 


———* = — pert Wr, (6.38) 


_ 3a Ude T, Om; 
rn = 539 OS (145) (aati) 9 


Thus, if one neglects the turbulent energy outside the energy-containing 
region connected with the maximum in the spectrum, one can put 


Q 
Ver 


where 


W= (6.40) 


The dependence of y.q on the energy input Q changes the picture of the 
turbulent spectrum, that is, increasing Q increases yj and it, therefore, 
can be shown that the maximum in the spectrum does not change when 
Q increases. The position of this maximum is determined by the relation 
T.\!4 1/2"-1) 
2a +7, Mee 
ol peg (ea cr (6.41) 
ks = 8m,c® : 


If the condition kKyjna,/k, < 1 is satisfied, as we have assumed, we have 
ko smaller, but not very much smaller, than k,. Thus, the position of the 
maximum is not very far from k,. 
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Another possible emission mechanism may be due to the processes 


I+iztti’, (6.42) 
or 
l+s = t. (6.43) 


The first process may be essential in an isothermal plasma when 
T, = T;, and the second one dominates over the first one in a non-iso- 
thermal plasma when T, > 7;. Both mechanisms give emission near the 
plasma frequency w,,, but the power emitted depends in a very essential 
way on the plasma inhomogeneity. In a homogeneous plasma both 
process (6.4) and the processes (6.42) and (6.43) transfer the energy 
of the l-waves as well as of the t-waves to small k-values, that is, to the 
region where the dielectric constant is zero and, therefore, such waves 
have difficulties in getting out of the plasma. For very small k the trans- 
verse waves have a frequency 

1 kc? 


Wo = Opty ——, (6.44) 
pe 


which is very close to w,., and then can be considered to be transverse 
plasmons in the sense that they are more like the eigen-oscillations of a 
plasma than like external electromagnetic waves. The optical depth for 
the processes (6.42) and (6.43) depends in an essential way on which 
part of the spectrum the waves come from which take part in the process. 
At the maximum value of k in the plateau region of the spectrum the 
optical depth of the processes (6.42) and (6.43) is of the same order of 
magnitude (we denoted it by Ly’**) 


max ~ © (Mi ae Opelole 
Lge = (7) a): (6.45) 
and at the minimum k-value on the plateau we have: 
min © (Me\™?  [ (@peroTe 
Lp x (7) 7). (6.46) 


Both of them can be less than the value L,, given by equation (6.33). 
If the plasma is optically thick for the processes (6.42) and (6.43), the 
turbulent energy on the plateau is of the same order of magnitude as in 
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the case when these processes are not present, but when there is a sub- 
sequent energy exchange between the longitudinal and transverse plas- 
mons. In the region of the maximum this kind of exchange does not change 
the total energy of longitudinal and transverse plasmons: 

W =W!'!+W* = constant. (6.47) 


6.6. Stochastic Plasma Heating in the Case of Langmuir Turbulence 


We established in the general theory that the stochastic heating of a 
plasma is due to the induced scattering of ions and electrons by the 
turbulent oscillations. The turbulent heating is due to the energy transfer 
in the region of damping. In the case of Langmuir turbulence, the energy 
is transferred from the region of Landau damping. Therefore, there is 
no turbulent heating in a homogeneous plasma. A change in the direction 
of the energy transfer is possible, if the plasma particles have an aniso- 
tropic distribution. This might be possible in magnetic mirror machines 
or in Q machines when the longitudinal plasma dimension is larger than 
the perpendicular dimension and the waves propagating in the perpen- 
dicular direction experience a density gradient which can change the 
direction of energy transfer. 

In an isotropic homogeneous plasma only stochastic heating is possible 
for Langmuir turbulence. Because the scattering by electrons is effective 
only in a very small k-interval, the ions are mostly heated stochastically. 
The largest growth-rate of heating is due to the largest k-value in the 
flat part of the turbulent spectrum. One finds that (Q,, is given below 
by (6.49)) 





oT; me\?5 Q T.\2 ] » if O<Qa; 
tf ~14(7*) : (147%) | (2 : 
i 9 1 ~ io 2) » if O>Qe. 
(6.48) 


The result (6.48) is valid for T, S 37;. In the opposite case, the factor 
Te]? /me\2!5 
[+a] Ga) 


is not present when Q > Q,,: the stochastic heating of ions is much 
larger than for the case when Q « Q,,; the quantity Q,, given by the 
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equation 
Ve m; \ ®/5 
oe — 6.49 
Ocr Ope MeNole ( Me 3 ( ) 


is rather small and usually QO > Q,. 


6.7. Stochastic Acceleration of Fast Particles 


Particles with velocities much larger than the thermal electron velocity 
can be resonant with the Langmuir oscillations, and therefore gain energy 
from them (Tsytovich, 1962). The non-linear effect on these particles is 
unimportant unless the resonance is allowed. This indeed becomes 
possible for waves with phase velocities larger than the velocity of light. 
Mostly they correspond to oscillations in the region of the maximum. Thus, 
generally speaking, resonant particle acceleration is due to the plateau 
part of the turbulent plasma spectrum. Sometimes in a sufficiently cold 
plasma (T, < 25eV for hydrogen) the tail of the spectral maximum can 
contribute to the particle accelerations. For non-relativistic particles with 
a kinetic energy e,;, the acceleration by the plateau region is described 
by the equation (€,;, = €—mc?) 





dexin _ Lr 
ieee ore 
where 
»_ 3 YBx) og mePmiPZ? () Te Q 
P= —"g "Peete ane (tn) Veo) (2) 


(m is the mass of the particle which is accelerated). In the region of the 
tail of the maximum, one finds an acceleration which increases rapidly 
with increasing energy. When the energy of the particle reaches the rest 
mass energy (that is, when the particle reaches relativistic energies), the 
diffusion coefficient in energy space becomes constant, and the accelera- 
tion follows from the equation e? = 4Df, or e/t = 4D/e; one finds that 
the energy growth-rate de/dt decreases with increasing energy as e~!. The 
characteristic e-dependence of de/dt is shown in Fig. 25. 

If the energy gain given by equation (6.50) is balanced by the collisions 
of the fast particles with the cold plasma, one finds for fast particles a 
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Sketch of the acceleration rate de/dt as a function of « for electrons and 
for ions in Langmuir turbulence (7 = 3v3,/cvy,). The part proportional 
to e%/4 is due to the spectrum k~4?, The part proportional ¢~-'/? is due to 
the plateau region of the spectrum, and the part proportional to e’-»? 
is due to the tail of the main maximum in the spectrum. If ¢ >> mc?, the 
acceleration rate does not depend on the particle mass and is proportional 
to e~1, The region de/dt ~ «("-"/? disappears when 7 > 1 (dashed curves 
a’ and b’) and the region de/dt ~ «7/4 when T, > T; (curves a and b). 


relativistic Maxwellian distribution with a temperature 


Ts = me y (a) (6.52) 


8? vn Te 
27x[1+(Te/Ti) ?memiWpe 


where 


On= (6.53) 


Usually Q>> Q,,, which means that 7, >> mc?. Thus the Langmuir 
TSY 7 
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turbulence acting for a long enough time creates relativistic particles, 
which in cosmic conditions can be identified as cosmic rays. In laboratory 
conditions the maximum particle energy is determined by the possibility 
of their confinement by magnetic fields, or, in other words, the maximum 
energy is of the order of the energy for which the gyro-radius of the 
accelerated particles is of the order of the characteristic installation dimen- 
sions. 

The number of particles accelerated is determined by the injection that 
can take place from the tails of the Maxwellian distribution. The most 
important role is played by low-frequency magnetic-type oscillations as 
injection mechanisms which can give preferential heavy-ion injection. 
When the number of accelerated particles is sufficiently high, the damping 
of the turbulent motion due to the acceleration can be important. This 
damping can exist approximately only on the plateau region, or more 
precisely for oscillations with phase velocities less than c. The turbulence 
can come to a stage in which practically the whole power is dissipated by 
fast particles. Under these conditions the maximum of the turbulent 
spectrum must disappear, because the turbulent energy must be absorbed 
by fast particles until it comes to the region where such absorption is not 
possible. The subsequent changes in the turbulent spectrum when the 
density of the accelerated particles is increased is shown in Fig. 26. 


W, 
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Change in the Langmuir spectrum as the density n4 of accelerated par- 
ticles increases; nt < n} < ni. 
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6.8. Correlation Effects for Langmuir Oscillations 


If there is a maximum in the spectrum, the correlation effects for 
Langmuir oscillations are mainly determined by the plasmon-plasmon 
scattering and the correlation width is given by the estimate 

1 wow, (W\? 
Aw ~ =~ Woe = [—-) - 6.5 
© * 30 P* v8, (r.) G2") 
This width is the half-width of the correlation curve near resonance. 


Very far from the resonance the correlation tail has the form shown in 
Fig. 27. 
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The correlations in Langmuir turbulence at low frequencies; 
Omar = 3k?V%./20 50. 


In this spectrum there is no relation between w and k, but the maximum 
frequency in the plateau region of this spectrum is equal to 


3k*v}. 


®max = Fy 
3a 
pe 


that is, the difference of two Langmuir frequencies. Such an effect was 
observed by Fainberg and coworkers (see Fainberg, 1967) in the beam- 
plasma interaction. 

The relation 





(6.55) 


A@max a ko 
aos aaa } 
®max k 





(6.56) 


and the existence of a maximum in the curve of the correlation function 
is due to the maximum in the turbulent spectrum. This maximum disap- 
pears when there is no maximum on the spectrum. 


q* 
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Results of a numerical computation of the correlation curve near w = w,. 

The different curves correspond to different k-values in the spectrum; 

k «k,; the parameters used in the computation are the same as for Fig. 

22. The non-linear shift in @ of the curves is small compared to the 
thermal one described by equation (6.1). 


In Fig. 28 we show the results of numerical computation of the correla- 
tion curve near w = c, for Langmuir turbulence. This figure represents 
the computation by Makhankov and Shchinov (1970) of the non-linear 
integral equation for the correlation function given above in Chapter 4. 


7. Electromagnetic Properties of a Turbulent Plasma 


THE presence of turbulent collisions has an important influence on the 
electromagnetic properties of the turbulent plasma response, when a 
weak field acts on the plasma. This means also that the growth-rate of 
several plasma instabilities can be changed by exciting turbulence in a 
plasma. This kind of change may lead to some stabilisation properties, 
as well as to new kinds of instabilities. One expects that this change in 
the electromagnetic properties can be essential for perturbations or for 
an external field with frequencies much less than the characteristic fre- 
quency of the turbulent collisions. This effect can be described by a 
dielectric constant depending on the turbulent energy. If the frequency 
is less than the turbulent collision frequency, one may expect that the 
dielectric constant cannot be expanded in terms of the turbulent energy. 

To show how such a theory can be developed, let us return to the 
general description of the turbulent state, introducing the stochastic and 
regular variables 


ae re Te es = / stoch, ana x 
-L— +(0-vf O+e(E B. le ~e( (B h ap )): (7.1) 
3) pigeon ogee a | te(Ee. =.) 


ar 
+e(Eses. a ae) -e¢ (Es =) =0. (7.2) 


Let us assume that in the stationary turbulent state the regular field 
E*® is absent. Denote all quantities for this state by an index “0”. Thus 


we have 
'stoch 
oH +e-vir) =—e (Beem EY, ay 
ofgion® toch Stoch , ofe** 
AE + (0-vfae)+e( Ege. 
+e( Egtoch, aorn ef Ejtorh. ae AG gs 0. (7.4) 
op \ / 
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Suppose now that there exists a small perturbation connected with the 
field E'*®, that is, 


fre — fice + Ares, fetoch = DP ieee + For, 
Estoch — Each + E3toch (7.5) 


where fyxet << fee; fStorh < pstoch; Estoch < Fstoch and ail quantities with 
index “|” are proportional to Es, 
From equations (7.1) and (7.2) we get the system of equations 


ofr® reg re, ofe"® 
Or +(0-Vff +e(E Be. op ) 


Soce (5 : CE pornyatonn + Ege fee) =, ; (7.6) 


stoch stoch ren 
ie Hove) be( Be. OS _Ofe" i ) +e(zamn. 2 ) 


+e Estoch, Of nay] eee LE sch psoch. 4. Batoch fatech 
1 op ~ 1 0 0 1 
—Cepstfeens Egeenfgo®)) = 0. (7.7) 


One must add to these equations the Maxwell equations (for the sake of 
simplicity we have assumed that the fields are longitudinal) 


div Ere? = 4 pee [amg “se yi (7.8) 
3 
div Ez = dn Ye, | fploch Spe (7.9) 


These equations are exact equations which describe the perturbations 
from the stationary state. Various approaches similar to the ones which 
we used earlier to describe the stationary state can be used here for the 
perturbations. 


ELECTROMAGNETIC PROPERTIES 95 


7.1. Expansion of Turbulent Collision Integrals in Terms of the Turbulent 
Energy 


The collisions with the turbulent oscillations are, for the perturbation, 
represented by the right-hand side of equation (7.6), which we denote 
by J. Let us use an expansion in the turbulent field E3'° and calculate 
the first term of this integral which is proportional to (ESE Ss) = 
1,6(x+'). Thus in this case we consider only the first-order turbulent 
collisions, which are proportional to the turbulent energy, that is, the 
first term in the expansion in terms of W/nT. To calculate this integral 
one needs only the first-order term, that is, the term which is proportional 
to the turbulent field in E$'°°, f$'°°", and fete". 

Firstly, from equation (4.10) it follows that f3'° can be written in the 


form 
toch __ e stoch , Ofo 
RP = ( an. 2), (7.10) 


where f;'0°" is the Fourier transform of f;'°", and x = (k, @). 
Secondly, neglecting the non-linearity in the stochastic field term in 
equation (7.7), one finds for the Fourier transform of f;'° the equation 


stoch e _ stoch , Ofo® a. aed — 
1 i[o—(k-v)] \ >” Op i[w —(k-v)] 





x (2. i | dn, dit 5(% —x1 — x2) {Ene foo + Boe fre, ) (7.11) 
On the other hand, from equation (7.9) we have 


Eggs = “te | peicen (7.12) 


We have, for the sake of simplicity, taken into account only one kind 
of charge (say, the electrons) assuming that the average charge density 
is compensated by a neutralising background (of ions) and that this 
constant component does not contribute to the higher Fourier components 
considered in equation (7.12). It is easy to take other charge species into 
account in a similar way. 
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Inserting expression (7.11) into equation (7.12) one finds that 


ae ane | ds, d'x2 6(% — 1 — Xe) oni 


Korte (GZ lest testes), 7.13) 
where «, is determined by equation (4.21). Two points need to be im- 
proved in this analysis. Firstly, one cannot in general divide by w—(k-v), 
if it is possible that w can be equal to (K-v). Secondly, one cannot divide 
by «¢,, if ¢, is close to zero. These are the same singularities as in the 
description of the usual turbulent state. 

Nevertheless, if one introduces expressions (7.10) and (7.11) into the 
right-hand side one finds that 


hex = -(5, [a's , d'x ” 8(%— ee ar EE ee Boo fee stoch 


=, 2 4° 4% = ae ” stoch _ 1 = 
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D7 yore! Oc! ofe* stoch r \ 
En) BE) oo (RECN) ges . (7.14 
| Winn ee ba 0, x2 op; +( 0, x1 Wh 2 / ( ) 
Bearing in mind that 
Stoc! SOC: k; ik; ” 
(Ba): Ease") = ~ pa be Oe" +0), (7.15) 
and defining J, by the equation 
(ess) (E82)) = FEE 1, (7.16) 


we can rewrite the collision integral (7.14) in the form 


afe® ie 5, oe" 
I, = Din 
: Lae ”* Op; + Esl Din - op 
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where 
(x —x' x") kik;’ kiki’ 4 
(1) = ie? aaa ay se DN 4! Thal 
Diy, x ie { kk" Ie, x Lom =k ~B) To —(K-v) am dx , 
(7.18) 
kiki; dx 
(0) only lil; 1 
Di = i | ras PTB C19) 
i ke" @+o,—(K+K,]-0) +15 Op, [o1—(k,-v)—i5] 
(7.20) 


To finish these calculations, it is necessary to find J{?,... Multiplying 
equation (7.13) by E3'0" and using equation (7.10), one finds that 


_ Ane*l ! 
teen een ” 
Ie ape (Qn = 6(x' — x +x ee ey (0) 


in’ Ex! 
0 ae (t we, k" | : 
x (2 ° ae k" Be ap + fim d m1. (7.21) 


Introducing expression (7.21) into equation (7.19), one finds that 














Df, = ys (EX), DR. + Dis. (7.22) 
where 
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98 THE THEORY OF PLASMA TURBULENCE 


Finally the equation which describes the perturbation (7.6) can be 
written in the form 


— ilo — (ken) fgtite (es. are | ES ae (7.25) 
whcre the collision integral has the form 


oft’ ) e 
El De Tt Bp De Gp, th OEE) 


a at 3} 
xi {Birt Din} - a (7.26) 

This equation is not as simple as it seems, because the D{P, contain 
the perturbation f{*. Therefore, equation (7.25) is an integro-differential 
equation for f°. The physical meaning of each term in equation (7.26) 
can be easily understood: the first term describes the change of the quasi- 
linear interaction, if such an interaction is possible. The other terms are 
negligible if the resonance w = (k+v) is possible. If @ 4 (k+v) the third 
term describes the change in the non-linear scattering and in the decay 
process. ; 

Because the w and & of the wave considered are now arbitrary, one 
can say that «,_,, is close to zero, that is, one can consider the frequencies 
which are close to differences in the frequencies of the turbulent oscilla- 
tions. This shows that the contribution from the non-linear and Compton 
scattering 1s small. Thus one must consider mainly the case of resonant 
turbulence with coefficients D®© and non-resonant turbulence with 
coefficients D™’, 

Before starting on a detailed analysis of equation (7.25) let us emphasise 
one important feature of this equation, that is, that even if we expand 
the J, in terms of the turbulent energy we find that the dielectric constant 
of a turbulent plasma cannot be expanded in terms of the turbulent 
energy. Indeed, if the frequency @ (or more correctly w—(k-v)) is less 
than the turbulent collision frequency, which can be estimated from the 
right-hand side of equation (7.25), the plasma response to the field is 
determined by the turbulent collisions which appear in the denominator 
in the dielectric constant. This also can be illustrated by the well-known 
dielectric constant due to the ordinary collisions in the region where the 
perturbing frequency is less than the collision frequency: 

wee 


ex l—-- . 
lOVeon 


(7.27) 
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Very roughly speaking, an equation like equation (7.27) occurs now 
where instead of ».,, an effective turbulent collision frequency appears, 
which depends on W, so that equation (7.25) is proportional to the 
turbulent energy. Thus, ¢ of a turbulent plasma cannot be expanded in 
terms of the turbulent energy. 


7.2. Effects of Turbulent Resonance Broadening 


The broadening due to the wave-particle interactions in /3'°" has 
already been discussed and thus instead of equation (7.10) one finds 


ofr 
Ste Para ieg(w = (k+v)) (zs < *) 7 (7.28) 


A similar technique can be applied to equation (7.7). Because the non- 
linear terms occur only in e([0/Op]-[Eg'* fo" — (Ege f°) ]})_ the 
interaction procedure leads to the same substitution of g for [w—(k-v)]"! 
in equation (7.13). Other interaction terms are of higher order in the 
stochastic fields. For small @ and k the factor e71,, in equation (7.23) 
becomes large because ¢,, ~ 0. One needs to take into account the non- 
linear corrections proportional to W,. This was done by the author 
(Tsytovich, 1969b). Roughly speaking, the result of such a consideration 
is to change «, to e, ten". 


7.3. The Dielectric Constant for Ion-sound Turbulence 


The most important interaction for the case of ion-sound turbulence 
comes from D{?, (see equation (7.19)). One can consider two limits. 
First, the case when w « w, and k « ky, that is, the perturbation with 
frequencies much less than the turbulent frequency, and wavelengths 
much longer than the turbulent wavelength. In that case, the diffusion 
coefficient, 





pw ane [ fs Bailey 1, 5(0—(ky+0)) d*xy, (7.29) 


does not depend on @ and k and is the usual quasi-linear diffusion coef- 
ficient. This means that at low frequencies the turbulent collisions act 
as real collisions. 
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Secondly, the case of high frequencies and short wavelengths, w > au, 
k > ky; in that case, the diffusion coefficient is proportional to 6(@—(k+v)) 
which means that, if for certain values of w and & the resonance condition 
is not fulfilled, that is, @ + (kv), as for example happens for the trans- 
verse waves, the diffusion coefficient is zero. (Of course for transverse 
waves it is necessary to generalise equation (7.19) to take into account 
the transverse field, but the conclusion is the same.) Thus the high- 
frequency waves do not “feel” these turbulent collisions. This shows the 
difference between the turbulent and ordinary collisions.* 

Equation (7.25) can be written for the case of small @ and k in the form 


ak ‘s reg reg, Ofo° ae 9) oft 
ileo—(k-v)) fi tte E: ap) = y 5, ae (730) 
The 6-function in equation (7.29) is only an approximation of the g-func- 
tion which broadened the resonance. 

The presence of the collision term in equation (7.30) rather changes 
the electromagnetic properties of a plasma in the low-frequency region. 
One can say the turbulent Ohmic heating changes the properties of the 
drift waves. If the spectrum of the turbulence depends only on |k,| and 
on x), which is the cosine of the angle k,; makes with some chosen 
direction (for example, the direction of the electric field), then the 
diffusion coefficients Df) can depend only on {|p| and on & which is the 
cosine of the angle p makes with the same direction. Because angular 
scattering is the most important effect for ion-sound turbulence, in the 
right-hand side of equation (7.30) only 0/0& appears: 


— ilw —(k+v)] fy°e+ eB : ) 


LD re? { x81, 5((Ky-0)) dae, et 31 
=> PF sme xj a (¢ 1*0)) d Hy OE . (7. ) 


Equation (7.31) can be rewritten as follows: 


eas re; re ofo® 
ifco—(k oMfise+e( Bs oe 


2 Vva-&) 2 
| x | _ dx i We, xo deo LE (732) 





PPVWpi V(1-&—x*) OE 
One finds thus that the collision integral is just the same as the one which 


t For high frequency waves there still exist anomalous collisions, but they are 
less important than the quasi-linear ones. 
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describes the collisions of the unperturbed electrons with the ion-sound 
waves. The difference lies in the possible polarisation of the field E;*® 
which can be non-longitudinal. Because fo is approximately isotropic 
we have 

ofge® pais ae 


ap = arg (7.33) 


If the turbulent collisions occur often enough, one can neglect w—(k-v) 
on the left-hand side, and then by averaging over the angles perpendicular 
to the z-axis one finds clearly for E/** the same anomalous conductivity 
as for the longitudinal field. Thus we have 





( iD), = | ev fics d3p = gtr Eres) (7.34) 
J zwi1,x* (2x8 * z? 


where a> is the turbulent conductivity. The other components are 
slightly different, but of the same order of magnitude. One obtains the 
following estimate for the skin-depth for an electromagnetic wave entering 
the turbulent plasma: 

c 


7.4. Electromagnetic Properties of Langmuir Turbulence 


Neglecting the renormalisation of the plasma Green function one 
finds the dielectric constant which describes the plasma response to the 
field EX*® as follows: 
w2. 


kvFel 1 + (r9/V5-)4) ” 


where w « kuz, © «K @,, k « k,, while ef? is a normal linear response. 
The coefficient d has the following form: 


e, = e+ (7.36) 


4x 0 
d=-, Tie i Pha Fay («. or ai) ¥> (7.37) 


where U,, = dw,/dk, is the group velocity of the Langmuir waves. The 
instability found from this dielectric constant corresponds to the result 
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first found by Rudakov and Vedenov (1964) and by Gailitis (1966). 
Such an instability can exist only if one can neglect the non-linear 
renormalisation of the plasma Green function. Analysis shows that it is 
usually difficult to satisfy the conditions for the neglect of this non-linear 
renormalisation. For example, the total turbulent energy must be con- 
centrated in a very narrow interval of the turbulent energy. For the distri- 
bution found above for the Langmuir turbulence spectrum, it is very 
difficult to satisfy these conditions in practice (Tsytovich, 1969b). There 
still exists the possibility to have a transverse perturbation described by 
the transverse dielectric constant ¢-. 

If |e'] >> 1, we have |H| ~ +/(\e'|) |E| >> |E| which means that the 
magnetic field is excited, if the plasma becomes unstable. 

This kind of instability is similar to that of magnetic field excitation 
by a conducting turbulent liquid. The growth-rate y for the excitation 
of a magnetic type perturbation by Langmuir turbulence was found to 
be (Tsytovich, 1968b): 


p= ie (z) WV +4 / (3) es) : i, (7.38) 


il t\ v w? 
amie pe Te 3 Pe 
Ux = 12 / (3) No = { Wi, d ky ke 2 . (7.39) 


Some inequalities must also be satisfied which give the range of para- 
meters for which this instability can be developed. 

One can find also the skin-depth for a low-frequency wave penetrating 
into a semi-infinite plasma which in the case of w « »—_ is described by 


where 








ce Ver 
aed ae (7.40) 
where 
_ Wi, d "kw. 
ra =o (| “tat: ) ay 


Usually the increase of 6 by turbulence is more efficient for ion-sound 
turbulence than for Langmuir turbulence. 

One can thus see that the turbulent collisions give plasma characteris- 
tics, which can easily be measured. 

The Langmuir turbulence oscillations were shown earlier to be con- 
verted into transverse waves with frequencies very close to the plasma 
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frequency. This looks like a transverse plasmon. The presence of a trans- 
verse plasmon can give rise to an instability which is similar to second 
sound (Tsytovich, 1969c): 

y = kun, (7.42) 


1 Ww 
pers / (sears): (7.43) 
e 


This type of instability can be excited also by external transverse waves, 
if their spectrum is sufficiently broad. Usually, some frequency broad- 
ening arises even if the initial signal is narrow in frequency. This can 
be the case when the wave propagates in the direction in which the mean 
plasma density is increasing until it reaches @,,. On the other hand, a 
high-frequency field excites Langmuir waves. When its frequency 
reaches the region where it can be equal to the Langmuir frequency, 
Wpet+ $k vi_/pe, the excitation starts from k close to Ap!, when the 
Landau damping becomes unimportant. This is the frequency which is 
approximately equal to (1+3)@p¢. Then, if the development of the 
instability is faster than the characteristic time for a frequency change 
due to the wave propagation in the direction of the density gradient, an 
appropriate amount of the wave energy can be transferred into Langmuir 
waves. This energy then comes back into transverse waves, and so on. 
This process gives a subsequent non-linear lowering of the electromagnetic 
wave frequency. In order that transverse plasmons are present, it is neces- 
sary that the plasma is optically thick. The presence of Langmuir waves 
and especially transverse plasmons can change the properties of drift waves 
(Krivorutskii, Makhankov, and Tsytovich, 1969). This can lead to some 
kind of drift-wave stabilisation. For example, one can find that kinetic 
drift instabilities can disappear when 


_ OlnT, 
"= Sinn 





> 2. (7.44) 


This type of stabilisation is different from the high-frequency stabilisation 
only in the sense that the high-frequency field has a stochastic nature, 
and, therefore, a broad spectrum, and also the distribution of this field 
is not arbitrary, but corresponds to the self-consistent picture of the 
plasma turbulent spectrum. 

It is also possible to stabilise some of the hydrodynamic drift-wave 
instabilities. 
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For ion-sound turbulence the presence of turbulent collisions plays a 
role similar to ordinary Coulomb collisions, as they have the same 
velocity dependence. But the corresponding %%, is, of course, very much 
larger then the Coulomb collision frequency and is also anisotropic, if 
the ion-sound turbulence is anisotropic (as in the case of a current-driven 
instability). As a rough approximation, one can estimate the drift 
instability growth-rate in the presence of ion-sound turbulence by 
putting % instead of »,,). Usually dissipative drift instabilities develop 
when ion-sound turbulence is present. 

We thus see that the presence of turbulence can have an essential 
influence on plasma confinement. 


8. The Cosmic-ray Spectrum 


THE plasma that can be found in cosmic conditions is often turbulent, 
and the turbulent plasma state may be assumed to be a natural one in 
the universe. The great efforts which are being applied to trying to 
stabilise plasmas only emphasise this statement. One can then ask the 
question: what kind of consequences can be deduced from this point 
of view? The most important of these seems to be a production (or, 
rather acceleration) of fast particles. From the theoretical treatment 
given above, it follows that the presence of fast particles in the turbulent 
regime is a consequence of the presence of turbulence. This means from 
the point of view of energetics that the creation of cosmic rays (or 
relativistic particles) by a turbulent plasma does not meet with any essential 
difficulties. It is also well known that in galaxies there are enough turbu- 
lent energetic sources which can give energy to the cosmic rays (see Ginz- 
burg and Syrovatskii, 1964). Such sources may be supernova explosions, 
active nuclei of galaxies, and so on. It is also known that cosmic rays 
are produced by solar flares, which represent an example of an explosion 
on a small scale. The problem of cosmic rays is thus not so much the 
problem of energy being transferred to fast particles, but the problem 
of the distribution of this energy between them, or the problem of the 
cosmic-ray spectrum. The spectrum of cosmic rays produced by solar 
flares is different from that produced by radio sources such as supernova 
remnants and radio galaxies. The solar flares have a very steep spectrum 
which can be represented by e~? where y(~ 7 to 8) is a large number, or 
by an exponential spectrum e—*/*o, The last spectrum exactly corresponds 
to the predictions from Langmuir turbulence, which creates a relativistic 
Maxwell distribution of fast particles. The radio sources have a power 
spectrum e~” where y usually lies in the range | < y < 3 with an average 
value of y of 2-5 to 2-6. These are the cosmic electrons. The cosmic 
ions—mostly protons—have a power-type spectrum near the earth, with 
y ~ 2-7. Thus the mean energy of the particles lies near mc?, but there are 
a large number of very energetic particles. The question is then how the 
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energy of the particles can be distributed in such a way that a power-type 
spectrum is created and why the spectrum of solar-flare cosmic rays is so 
different from the spectrum of radio sources. One may assume that the 
difference is connected with the dimensions of the sources, and we shall 
see later that it could be connected with the optical thickness for the 
plasma emission mechanisms. 


8.1. Energy-dependence of the Acceleration Rate 


Let us consider the energy dependence of the acceleration rate. Suppose 
for simplicity that 


de 
“ai — Ben “ (8.1) 
The three different cases, « > 1, w = 1, and ze < 1, are shown in Fig. 29. 
de 
dt p.>1 


w= 


<i 


€ 
Fic. 29. 
The acceleration rate as function of the particle energy. 


The characteristic acceleration time t can be defined as 


E 1 

t= deja = B el-#, (8.2) 

If « < 1, the larger the energy, the longer the time needed for accelera- 
tion. Thus, a particle injected with ¢,;, into the acceleration regime and 
reaching &max needs to be accelerated in a time that is determined by ena. 
The opposite occurs for x > 1. In this case, the larger the energy, the 
shorter the time needed for acceleration. This means that the maximum 
possible energy is reached in a time determined by the injection energy. 
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The case gs = 1—which, for example, corresponds to the Fermi 
acceleration—is intermediate. In this case, the time is independent of the 
energy (more precisely, logarithmically dependent). 

The difference between the case : < | and uw > 1 plays an important 
role in the formation of the particle spectrum. When « < 1 the particle 
spectrum, that is, the dependence of the particle distribution function 
f, on the energy ¢, must have a very important cut-off at sufficiently high 
energies. For example, in the case of the Langmuir turbulence considered 
above, the law for the acceleration of relativistic particles is de/dt = «/e, 
and if this acceleration is compensated by ionisation losses, the distri- 
bution function f, has an exponential cut-off. Indeed, as was shown 
above, the distribution of fast particles forms a relativistic Maxwell 
distribution with a temperature 


Te >> me*. (8.3) 


This 7. plays the role of a cut-off factor in the distribution function. 
We shall normalise the distribution function of the relativistic particles 
as follows: 


{ "fds =n, (8.4) 
0 


where n is the total density of fast particles. Because these particles have a 
maximum energy much higher than mec? the non-relativistic part of the 
distribution function is sometimes of no interest. In this case, « = cp 
and the phase volume p? dp is proportional to ¢2 de. Thus, the distribution 
function with T = 7T,, can be written in the form: 


fe po ere Tet ¥ (8.5) 


This distribution function seems to be different from that which 
corresponds to the radio-emission in radio sources. This is the case, 
because the power-type emission J, ~ w~”" seems to be possible only for 
a power-type particle distribution. To form a power-type distribution 
from the distribution (8.5) is difficult. Indeed, the conditions in different 
parts of the source must be different. The radiation observed is the super- 
position of the radiation arising in different conditions. Thus, one can 
average expression (8.5) over a distribution of T._-values. In order to 
have a power-type distribution one needs to have a power-type distri- 
bution for T.g which seems to be quite an artificial assumption. 


8* 
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This is one example of the consequences of the case when uw < I. In 
other examples which can be worked out, the distribution function cannot 
be exactly of the type (8.5), but it is difficult to believe that in every case 
with a strong cut-off the distribution of these will be such as to give a 
power-type spectrum. 

The sharp cut-off is perhaps consistent with the cosmic-ray spectrum 
produced by solar flares, but not with the cosmic electron and cosmic 
ion spectrum. It is also necessary to emphasise that in the case of Lang- 
muir turbulence the spectrum (8.5) arises only in rather dense plasmas 
with rather small values of Tog. As we shall see later, Langmuir turbu- 
lence creates a power-type spectrum for sufficiently high particle energies. 

In the case when pw = 1, the well-known example of the Fermi accelera- 
tion leads to similar difficulties. Indeed, it is known that if the time of 
acceleration is t, Fermi acceleration leads to a power-type spectrum 


with 


y=lt ae (8.6) 
where 
a= : - (8.7) 
Both the parameters « and t can vary in broad intervals so that 
l<y<o. (8.8) 
This contradicts observations which give approximately 
l<y <3. (8.9) 


This fact is normally the basis for rejecting the use of the Fermi accelera- 
tion to explain the cosmic-ray spectrum. We want to emphasise here 
that the greatest difficulty 1s that the average spectrum will not be a 
power-type one at all. Indeed, the real spectrum is due to the super- 
position of different types of spectrum produced in different regions with 
different « and t. Thus 





E 


f= | const F(y) dy, (8.10) 


where F(y) is determined by the distribution of « and t. Only for very 
special choices of F(y) will f, be a power-type function. Thus the Fermi 
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acceleration, as well as any other process with « = 1, meets a difficulty 
similar to the uw < 1 case. 

One of the reasons that y varies in such a broad interval as given by 
equation (8.8) is that the acceleration and energy-loss mechanisms are 
quite different. 

Thus we see that neither « < 1 nor u = | are appropriate for explaining 
the power-type spectrum. Thus only the 4 > | case seems to be possible. 
Of course, 4 < | is very effective for small energies, and can, therefore, 
be considered as a good mechanism for the injection. We shall see that 
the case « > | can start only at rather high energies, « > ¢,, where ¢, can 
be called the injection energy. If this happens the lowest acceleration rate 
is for ¢ ~ e€,. One of the interesting points of the acceleration for the 
case when 4: > | and ¢ > «, is that the highest energies of the cosmic 
rays are created very rapidly, and one can expect even now to have an 
acceleration of particles to the highest energies. This question is connected 
with the problem of forming the highest energies in the cosmic spectrum, 
which are of the order of 10?°eV. One can estimate that if the particles 
are protons the friction produced by the 3°K black-body radiation 
essentially reduces their energy on a time scale of 10® years. The question 
is then whether such particles can be accelerated in some regions on a 
very short time-scale. To have an answer to this question, one does not 
need a precise model of the source of the cosmic rays. If 4 > 1, one need 
only assume that cosmic-ray creation exists at all, that is, that the rate 
of acceleration is sufficiently high to compete with the energy losses at e,. 
If this is the case, one can also ask what is the highest possible energy 
€s» for which yg remains larger than unity. One can ask what the plasma 
density, temperature, the dimensions of the source, and so on, will be, 
such that «,, is very high. After this one can ask whether in cosmic 
conditions these requirements can be met near very energetic sources 
such as pulsars, quasars or galactic nuclei. Thus the problem of the cosmic 
ray spectrum becomes of a kind such that it can in principle be answered 
simply by theoretical considerations. 


8.2. Energy-dependence of Resonant Wave-particle Interactions 
The rapid rise with energy of the acceleration rate can only come from 


the rapid rise of the particle-energy diffusion coefficient. This takes place 
if the larger the particle energy, the larger the spectral interval of turbulent 
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where o denotes the turbulent wave. In the high-frequency region, one 
can write k = w/c, and 


fia: w° —(k°sv) ee w° —(k°s0) 
~ 1—(v/c) cos 6 me 1=(v]/0) 
= _ 22 7 — (Kk? ) if 2, 8 17) 
= nel wo ov), if e>>mce?. (8. 


There are also two possibilities as in the case considered above, namely 
w’ >> k°c or w’ « k°c. The first case means that the phase velocity of 
the turbulent oscillations is larger than c, while it is less than c in the 
second case. All turbulent oscillations, except Langmuir waves, have a 
phase velocity less than c. 

In this case, only the first-order resonance wave-particle interaction 
is of interest (see equation (8.13)), as it can be proved that the second- 
order effect is always small, if the broadening of the resonance is taken 
into account. 

On the other hand, it was shown that the energy of the Langmuir 
oscillations is mainly concentrated in the region where the phase velocity 
is larger than ec. 

The condition (8.17) then has the form: 


E 2 
ies 2( snc wo. (8.18) 


This means that the higher the particle energy, the higher the frequency 
of the transverse waves that can interact with the particles. In order to 
have uw > 1, it is necessary that the energy of the transverse waves 
increases rapidly with frequency. We shall see that this naturally arises 
because of a re-absorption of high-frequency radiation by fast particles. 


8.3. Acceleration by Low-frequency Turbulence 


Let us start by analysing the lowest possible frequencies of turbulent 
oscillations, which are due to magnetohydrodynamic and Alfvén oscilla- 
tions. Most of the energy of such oscillations is contained in the cold 
plasma, but not in the electric field. This leads to the conclusion that 
only magnetic types of turbulent motion can give an effective acceleration. 
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Indeed, the diffusion coefficient for the sound wave w ~ kv, cos 0 (wa <«< 
@,;) is approximately constant, if the particle velocity is near to c. 
Therefore, <2 = Dt = const ¢ or e/t = const/e, and « =—1. Thus, the 
most interesting case is va >> v,. The two types of turbulent waves (Alfvén 
waves, w ~ kv, cos 0, and fast magnetohydrodynamic waves, w ~ kv) 
accelerate in the same way. The distribution of these waves in & or in 
@ = kvy is shown in Fig. 30. 


Wa 





Wo = koa Oy ©) 
Fic. 30. 


The spectrum of collisionless magnetohydrodynamic turbulence; ky ~ 

1/L), where L, corresponds to the characteristic scale in the energy- 

containing region. Here kyv, is of the order the ion-ion collision fre- 
quency »;. 


The value ko corresponds to the energy-containing region. The drop 
near Wy; is due to the ion cyclotron-resonance absorption. In the collision- 
less region the energy flows in the direction of lower frequencies. The 
value ko may correspond to 1/Z, where Lj is a dimension of the plasma, 
or to @) = 1, where », is the ion-ion binary collision frequency, provided 
v,/Lh <« ¥,. In the region w « », the ordinary collision-dominated mag- 
netohydrodynamics lead to energy transfer to the highest frequencies. 
Therefore, one can expect for sufficiently large astrophysical objects that 
ky = v/v,. If vg is of the order of the ion thermal velocity, ko is of 
the order of the mean free path of cold ions in the plasma. The accelera- 
tion of particles due to the spectrum of Fig. 30 is shown in Fig. 31. Here 
Ex» 1S of the order of the energy for which the gyro-radius of the accel- 
erated particle reaches Lo; we find thus from equation (8.15): 

2, = ae VA. (8.19) 


v; 
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Fic. 31. 


The acceleration rate of particles which are accelerated by collisionless 
magnetohydrodynamic turbulence. 


The rise of the acceleration rate with increasing particle energy is due 
to the fact that the resonant condition (8.15) raises the number of waves 
interacting with the particles. The decrease of de/dt when é > ¢,, is 
due to the fact that all the waves in the spectrum are already included in 
the acceleration of the particles. The coefficient « is of the order of 


vA WwW 


“Loc H/8x’ 


(8.20) 


where W is the turbulent energy and is of the order of (H2 /8z) where H. 
is the turbulent magnetic field associated with the Alfvén and magneto- 
hydrodynamic waves. 
One can compare equation (8.20) with the Fermi acceleration coef- 
ficient: 
ap = u2/Le, (8.21) 


where uw is the velocity of the magnetic clouds and L the distance between 
them. One can see that v, plays the role of u, and L, the role of L in this 
comparison, and there is an extra factor (H?)/H? which is usually « 1. 

The acceleration by magnetohydrodynamic and Alfvén waves can, 
therefore, be considered equivalent to the Fermi acceleration in modern 
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turbulent plasma theory. One concludes then that this type of acceleration 
is always weaker than the Fermi type. This is so because of the factor 
(H2)/H? and also the factor (e¢/e,)*-”, if » < 2. If » > 2, we have w < 1 
and this does not seem to be a case of interest. The case » < | can be 
shown to be unstable for the formation of a turbulent spectrum (Lifshitz 
and Tsytovich, 1970). If 1 < » < 2, we have 2> p> 1. If w > 1, that 
is, » < 2, the factor (e/e,)?~” is always small, except when ¢ ~ e,, and 
the acceleration rate is less than for the Fermi case. The essential point 
is that Ly is very large, of the order of v,/v,, and, therefore, the « of 
equation (8.20) is rather small. There exists some speculation in the 
literature about the possibility of raising « in collisionless regions by 
decreasing Lo. In all these estimates, it is necessary to consider more 
carefully the possible time that such a spectrum with small Lo can exist. 
This time can be determined by: (i) non-linear spreading of the turbulent 
energy which raises Lo, (ii) absorption of the turbulent energy by cold 
plasma particles, (iii) absorption due to fast particle acceleration. The last 
possibility was considered by Tverskoi (1968) and can only occur when 
the first two possibilities do not work. One can estimate the characteristic 
time of the non-linear interaction to be of the order of 


1 ~ 
yal A Se kva . (8.22) 


If one substitutes here for k the minimum value which corresponds to 
equation (8.15), one finds 


1 vio me H2 
gi On ge OE EE (823) 


The non-linear interaction which raises Lo is unimportant only if the 
characteristic time of the non-linear interaction is greater than the 
acceleration time. Overestimating the last one as I/a« (for ¢ = &,) we 
find the following inequality, which shows when one can neglect the 
non-linear interactions: 


(ie eee (8.24) 


However, the minimum Lp is c/@,; (@ ~ w,;)and this means €/mjc? > 1. 
The last condition can be satisfied only for ultra-relativistic or heavy ions. 


THE COSMIC-RAY SPECTRUM 115 


On the other hand, « of the accelerated particles cannot be very high, 
because their Larmor radius must not exceed Lo if they arc to be efficiently 
accelerated. This gives: 





— boo Hi_ (8.25) 


at A 1, (8.26) 


This last relation needs extremely high magnetic fields. This means that 
one can find a turbulent spectrum where Lo is determined by the fast 
particle acceleration, only in the extremely peculiar cases of very high 
magnetic fields, and high particle energies. Even under these conditions, 
efficient acceleration can exist only in a very narrow particle energy 
interval, and the spectrum produced is not a power-type spectrum. 
Therefore, one cannot apply the results found from the theory of a 
spectrum determined by fast-particle acceleration to explain both the 
acceleration of non-relativistic particles in the interplanetary space and 
the cosmic particles in the galaxy. 

To be certain that the other processes cannot diminish Lo, one needs 
also to estimate the damping of the oscillations by the cold particles. 
This occurs mainly for the magnetohydrodynamic oscillations. However, 
the non-linear interactions very rapidly exchange the energy between 
the magnetohydrodynamic and the Alfvén oscillations. Both the damping 
and the non-linear interactions are proportional to w and this means 
that if the non-linear transfer exceeds the damping for one a, it exceeds 
it also for others. Thus the damping cannot form a cut-off of the turbulent 
spectrum at a particular value of w or be responsible for shortening Lo. 
The inequality which shows when the non-linear interactions are stronger 
than the damping has approximately the following form: 


2 2 
a = y me Us (8.27) 


m, VA 





If equation (8.17) is satisfied one can compare the fast-particle accelera- 
tion rate with the non-linear interactions and this gives relation (8.24). 
If equation (8.27) is not satisfied, the damping exceeds the non-linear 
interactions and the turbulent spectrum is not formed. One can, of 
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course, have an instability in a broad frequency interval, but then the 
energy of the turbulent oscillations is increased to the value for which 
the non-linear interactions become important. 

These estimates show that the spectrum of turbulence is mainly 
determined by the non-linear interactions and not by the particle accelera- 
tion. This means that Lo is rather large and excludes any speculations 
about diminishing Lo in the collisionless region. This means that we 
return to the results found in the paper by Ginzburg, Pikelner, and 
Shklovskii (1955) who estimated the acceleration of particles in supernova 
remnants; in this paper the Lo in the Fermi acceleration was taken as 
the mean free path of the cold particles (v, ~ v7; in this case). There is no 
possibility of applying the results of low-frequency magnetohydrodynamic 
turbulence to cosmic-ray acceleration in the interstellar medium. Thus, 
from the modern theory of plasma turbulence it looks as if Fermi-type 
acceleration mechanisms are rather unlikely. The conclusion that this 
acceleration is impossible in the interstellar medium, which one can find 
in the book by Ginzburg and Syrovatskii (1964), is even more conclusive. 

However, these types of turbulence, regulated by non-linearity, can 
give an efficient injection especially of heavy ions (that is, they give pref- 
erential ion acceleration) as was shown by Melrose (1969). 

The spectrum of the accelerated particles can be shown never to be 
a power-type spectrum, except in the peculiar case where v = 2. In this 
case, the acceleration is of the Fermi type and all the difficulties connected 
with it were already discussed. 

The reason why the spectrum is not a power-type one can be under- 
stood if we mention that the acceleration is due to the induced emission 
and absorption of magnetohydrodynamic waves, as the spontaneous 
emission is always negligible compared to the fast particle ionisation 
energy losses. Thus the acceleration and energy losses are due to quite 
different physical processes. 

If the acceleration by magnetohydrodynamic motion is not very 
efficient, one can ask what happens if one considers the different branches, 
or raises an effective Lo by considering the turbulent excitations of higher 
frequencies. The first of them is the whistler mode. Indeed such motions 
are known to be excited in the magnetosphere of the Earth, or in mirror 
magnetic confinement installations. One can find that at least these 
motions can give effective accelerations of non-relativistic particles, and 
for the ultra-relativistic particles the diffusion coefficient is approximately 
constant (Tsytovich and Chikhachev, 1969) and therefore up = —1. 
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8.4. Acceleration by High-frequency Turbulent Oscillations 


As was already mentioned, even the vj, <c part of the Langmuir 
turbulence spectrum gives an efficient acceleration of particles. The 
effective temperature of the particles easily reaches ultra-relativistic 
values. However, the spectrum of the accelerated particles is not a power- 
type one for e « Torr: 


S-® const -e?; exp (—e/Ter) = 1. (8.28) 


The normalisation of expression (8.28) is 
| fide =n. (8.29) 


The cut-off of the spectrum at 7. is due to ionisation losses. There exist 
conditions for which the change of the spectrum of the fast particles 
arises at energies much less than 7.,, so that the ionisation losses do 
not become essential at all. We shall see that this change leads to a power- 
type spectrum. 

First of all, we mention that the high acceleration rate by Langmuir 
turbulence oscillations is due to the fact that their frequency is rather 
high. We have already shown that the most effective acceleration by 
magnetohydrodynamic oscillations occurs at a frequency @ of the order 
of »;.On the other hand, the plasma frequency w,, is larger than that by 
at least a factor +/(m;/m,) Np, where Np is the number of particles in the 
Debye sphere, which in astrophysical conditions is very high. An idea 
of the order of magnitude of Np is given by the values in Table 8.1. 
The Fermi acceleration is due to the reflection of the particles from the 
magnetic clouds. The acceleration by magnetohydrodynamic waves is due 
to the change of the particle energy in the case when the particle is not 
reflected, but propagated through a layer with an increasing magnetic 
field. This effect is smaller than that due to reflection by a factor H2,/H?. 
The less the distances between these layers, the more frequent are the 
“collisions” of particles with the magnetic waves and the larger the 
acceleration rate. The same is true for Langmuir waves, which can be 
considered as consecutive charged layers with opposite signs of the charge. 
As in the case of magnetohydrodynamic turbulence, the non-linear effects 
increase the wavelength and decrease the acceleration. There exists, 
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TABLE 8.1. NUMBER OF PARTICLES IN A DEBYE SPHERE FOR VARIOUS PHYSICAL SYSTEMS 














Type of plasma nin cm-? T, in eV | Ny 
Thermonuclear 10%4 5x 104 10° 
Tokomak 2x 1018 1500 1-3 10? 
Laser 102 104 104 
Stellarator 1038 100 3x 10° 
Radiation belt 10? 4 2:6X 108 
Solar corona 108 100 108 
Interstellar medium 0:03 1 5x 101° 
Photosphere 105 0-5 10 
Cosmic rays 10-7 109 1076 
Pulsars 101° 109 3x 107 





however, an essential difference between Langmuir and magneto- 
hydrodynamic turbulence in that this acceleration does not decrease 
indefinitely. This is the result of increasing the phase-velocity of the 
waves which becomes greater than the light velocity. If this happens, the 
turbulent oscillations become non-resonant with the particles or, more 
exactly, Cherenkov resonance is not possible for them. In the region 
Vpn => ¢, Where most of the turbulent energy is concentrated, another 
process of scattering which is described by equation (8.17) becomes 
important. For the resonance Cherenkov interaction the minimum wave- 
length is thus the order of 


c 
Lnin © — - ’ 8.30 
pe — 
compared to 
Va VA m ¢ 
Dg es ee ee - 
aes - a Gia Np (8.31) 


for the magnetohydrodynamic turbulence. The increase of the power 
input Q in the Langmuir spectrum raises the Lo of this spectrum that 
determines the position of the maximum ko ~ 1/Lo (sometimes ko is not 
changed, as was already noticed). This does not affect the resonance 
Cherenkov interaction, though, because in most cases the phase-velocity 
of Langmuir waves which corresponds to ko is much higher than the 
velocity of light. 
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One must mention also that if the plasma has very large dimensions 
so that it is optically thick for the non-linear transformation of Langmuir 
waves into transverse waves, the conversion of longitudinal into transverse 
waves becomes important. This mainly happens when the frequencies of 
the transverse and of the longitudinal waves are equal, that is, wt = w!, or 


kc? 3 k*v3 
pet Dope = Ope + oY a ’ (8.32) 
le 


so that the frequency of the transverse waves are very close to the plasma 
frequency. In this process of conversion the frequency is slightly lowered. 
The energy converted into transverse waves is again converted into 
longitudinal waves. This subsequent oscillation of the energy between 
transverse and longitudinal oscillations for vp,,<e¢ was investigated 
numerically by Kaplan and Tsytovich (1967). An essential point is that 
these transverse motions usually have frequencies which are lower than 
those of the longitudinal modes, and they thus take part in the process 
of the frequency lowering. Therefore, they are so closely connected with 
the plasma properties that one needs to consider them as special trans- 
verse plasmons which are present in the plasma. The Compton scattering 
by longitudinal and by transverse plasmons does not depend on their 
distribution in k-space, as is obvious from the fact that the resonance 
condition (8.18) is independent of k for k<,,/c. Moreover, the prob- 
abilities for both Compton effects are equal and the result is independent 
of the subsequent transformation of the energy from longitudinal to 
transverse plasmons and back, in the region vp, > c, and depends on 
the sum of their energies, 

W=W'+W,—7 (8.33) 


only, and this sum is constant. The relativistic particles suffering Compton 
scattering while condition (8.33) holds are limited to high frequencies 
@ >> Wye as follows from relation (8.17). Their emissivity was found for 
power-type particle spectra, f, = const e~”, in a paper by Tsytovich and 
Chikhachev (1970). If Q,, is the spectral density of the power emitted, Q, 


Q= i Q., deo, (8.34) 
and if 


&, E 


7 taal (i) (8.35) 
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we have 
. a) (y+ D (y+ 3y+5) 
(8.36) 








This emissivity has the same dependence on the frequency as in the case 
of synchrotron radiation. Both effects are similar—synchrotron radiation 
is due to the curvature of the particle in an external magnetic field and 
this gives an emission with a frequency of the order of 


eH / « \? 
0 =< (raz) (8.37) 
while the plasma mechanism gives an oscillation with @,,. and emission at 
3 38 
ae (cr) (8.38) 


In the general case, both mechanisms must be taken into account. The 
curvature of the particle is not essential in the plasma mechanism if the 
wavelength is less than the gyro-radius, or 


ko> . ke (8.39) 


Only at this point does the Lo of Langmuir turbulence appear. The 
condition (8.39) can be written in the form 


ae Ope Me 
me ckyo m—~ (40) 
eH 1 
e pe 


As we shall see, this condition is usually fulfilled if £ < 1, as is normally 
the case under astrophysical conditions. This results also from the fact 
that a power-type spectrum is formed at rather high energies. 

The next problem is that of the re-absorption coefficient for this high- 
frequency radiation which is emitted. The emission considered above is 
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due to spontaneous emission. The balance of the induced emission and 
absorption leads, for an isotropic particle distribution, to a positive 
re-absorption coefficient: 


y ee Ope MET Oy 410) (Sy. @42) 
me nge,(y + 4) (y +6) . ee 

If the synchrotron emission is also taken into account, one must add 

the two re-absorption coefficients which are due to the plasma and to 

the synchrotron mechanism. If W is of the order of H?/8a the plasma 

mechanism gives much larger re-absorption, if € < 1. The equation of 

radiation transfer has in stationary conditions the form: 


I, 
+e(n-V71,) = O..-Yole s (8.43) 


where J, is the spectral density of the radiation. Therefore, we find (see 
Tsytovich and Chikhachev, 1969b) as a stationary and homogeneous 
solution 

1, = aw? , (8.44) 
where 


H?2 
se WOP+4y + 11) + EO aly) 


4 1 
a = Mbe ~ 5 (35 | ee es (BAS) 
1 Wt 6y+16)+ 60-2" arly) 2 





This kind of spectrum in the re-absorption region is well known for 
synchrotron mechanisms of emission and arises when both the plasma and 
synchrotron mechanisms are taken into account. The meaning of this 
is very simple. One can write J,,~ w?7eq; Te ~ € ~ </. Thus intensive 
transverse waves with frequencies much higher than w,, are created. 
As we have seen, the shorter the wavelength, the larger the acceleration 
rate. One can then consider the transverse-waves acceleration. The 
wavelength of the transverse wave is 


A= clo. (8.46) 


It is obvious that the higher the frequency the smaller A, but the spectrum 
(8.44) cannot be true up to infinitely high w because the total energy of 
the transverse waves then diverges. There are three possibilities for 
changing equation (8.44). Firstly, relation (8.44) is not satisfied for 


TSY 9 
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@ > ,, if the region becomes optically thin. If Z is a characteristic 
dimension of the system, this w, can be found from equation (8.42): 


clyoe = L. (8.47) 


The scattering of the transverse waves by the turbulent waves decreases the 
optical depth and leads to the system being optically thick for higher 
frequencies. The y,,, of equation (8.47) must be replaced by (see Kaplan 
and Tsytovich, 1967) 


yl = V/(YurFu)s (8.48) 


for the case when o,,, > y,,,- Here a,,, is a scattering cross-section which 
is approximately constant up to @ ~ (c/vze)@,- and which is of the order 


of 
3 
Sp) ee a 
o=4 / mt a Op / ( fat ) (8.49) 


while it decreases with increasing frequency as w~4, when w > (c/V7z_)Wre. 
The second possibility for changing equation (8.44) at high frequencies 
lies in quantum corrections. These are important when fw is of the 
order of e, or when 


2 
© ~ One . » if si ae (8.50) 


me me? —hpe 
or 
me ) 
O > Wg = Wpe{ ;---) - (8.51) 
a pe 
= / 


The quantum corrections are essential only if a, < ws, which can happen 
only in a very dense, large-sized and turbulent region. Usually the restric- 
tions on equation (8.44) are under astrophysical conditions due to the 
effect of the optical thickness; however, one cannot exclude the other 
possibility in the central core of galaxies or quasars. 

The third possibility is the non-linear interaction of transverse waves 
with one another which may become important as the spectral density 
I,, increases with frequency. This non-linear interaction has the form 

ol. ws oo, 


4! ayn. = 77 -__ PF at. 
a oe nm I 0m w- (8.52) 
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Introducing into y*" the J,, from equation (8.44), we have 





ws, ( w \1? 
VB! ~ py SP ta. ) (8.53) 


Comparing this relation with equation (8.42), we find the critical a = w"! 
for which expression (8.53) dominates: 


(8.54) 


al. E W ec i ( me? ‘a 3) 
ar" & Wpe ; 


No Mole UTe € 


* 


This quantity w™" is rather large, because of the factor Np which is 
usually very large. When w > w™! we have I, ~ @, and thercis no longer 
a rapid acceleration. 

The very rapid increase of energy of the transverse waves with frequency 
given by equation (8.44) in regions where it is valid leads to a very rapid 
increase of the acceleration rate with particle energy. Indeed, correspond- 
ing to relation (8.18) the higher the particle energy, the higher the fre- 
quency interval of the spectrum (8.44) which is involved in the accelera- 
tion process. This acceleration is due to the gain of energy by particles 
when the waves are re-absorbed. 

The spontaneous process and induced process of scattering of particles 
by the turbulence W must be taken into account; this was done already 
for the emitted radiation in equation (8.43). The same effect changes the 
distribution function of the particles: 


Be 4 (sift) = = ae De? . z +. - Af = 0. (8.55) 


The second term on the right-hand side of equation (8.52) corresponds 
to the spontaneous emission. One can calculate the energy change of the 
particles due to this term: 


Be fe © de eae tede | Anac 


ob i de | fae 7 fra 


Thus we have 


=A (8.56) 








ees, Be W. (8.57) 
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These losses have the same energy-dependence as the synchrotron losses 
and must be considered together with them: 


(8.58) 


2 2 
We ethes a noe e é ( A ) 


Ope et Eo ee. 
This term is balanced by spontaneous emission, that is, by the first term 
on the right-hand side of equation (8.43): 


At | OU: (8.59) 


The induced processes which in equation (8.43) are represented by the 
second term are balanced by the diffusion term in equation (8.52). If one 
calculates the quantity D in equation (8.52), one finds it proportional 
to 3, which results from the rapid increase in frequency of J,, given by 


equation (8.44). The quantity D is proportional to aa T,, dw/@, which 


is proportional to ¢°, if one takes into account that @ ~ «2. This e-depend- 
ence of D leads to an acceleration rate which is proportional to «2: 


“ = pe. (8.60) 


The same type of acceleration arises for synchrotron re-absorption, which 
was not considered properly earlier (more precisely speaking, the influ- 
ence of the re-absorption on the particle distribution). Both the synchro- 
tron and the plasma effects must be taken into account in equation (8.60). 
The stationary solution of equation (8.52) for the particle distribution 
function gives a spectrum of particles. If D ~ «3and A ~ e&2, this spectrum 
has the form e~”. Because all the coefficients included depend on y+ (see, 
for example, equation (8.45)) » must satisfy an equation of the form 


F(é, y, x) = 0, (8.61) 


where x = —In(W/nmc?). The solution of this gives y as a function of 
— and y. The result is shown in Fig. 32. 

All values of y lie in the interval 0:9 < y < 3 which corresponds to the 
observed values in radio-sources. There exists a special value of y which 
is of the order of 2-6 to 2-7, for which the change of W over 2 or 3 orders 
of magnitude does not appreciably change y». 
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Fic. 32. 


Results of a numerical solution giving y as a function of x and é. 


The reason why y is found to lie in a small interval is that both the 
energy losses and the particle acceleration are due to the same (sponta- 
neous and induced) process. One can expect under these conditions that 
the super-position of the spectra from different regions leaves the power- 
type character of the spectrum largely unchanged. In these calculations, 
the results were averaged over the directions of the magnetic field, 
assuming that it changes over distances which are large compared to the 
gyro-radius. The change of y from 0-9 to 3 corresponds to a greater and 
greater domination of the plasma mechanisms. Nevertheless, the y 
produced by the mechanism considered lie in the range corresponding 
to observation. 

Many problems must still be solved in order that the results can be 
applied to constructing radio-source models. The first of these is to take 
into account the particle anisotropy and to find a self-consistent solution 
in the presence of regular magnetic fields. This may also give a maser 
effect and polarisation of the emitted radiation. 

The second problem is to take into account the resonance diffusion 
coefficient due to particle acceleration by waves with w/k < c. This gives 
the connection of the distribution function (8.28) with a power-type 
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distribution function. The energy corresponding to this junction can be 
called an injection energy ¢, in this acceleration mechanism. An estimate 
for es can be found from the intersection of the curves for the acceleration 
rates de/dt due to resonance (w/k<«<c) and due to non-resonance 
plasmons, as shown in Fig. 33. 





Fic. 33. 


Acceleration rates caused by resonant (~ ¢~!) and by non-resonant (~ €?) 
plasma oscillations. The intersection point ¢ = €, corresponds to the in- 
jection energy into the fast-acceleration region (~ €?). 


Usually e+ >> mc? so that condition (8.40) is satisfied. The power-type 
spectrum disappears if es > Tyg, or if & > min {&s«, €4, e™!}, as we 
saw earlier. Thus one can find conditions for which a particle spectrum 
with a fast cut-off can be produced. In any case, to produce a power-type 
spectrum it seems necessary to have the system optically thick for relativ- 
istic-particle emission. This is not the case for solar flares and this may 
lead to a possibility for an explanation why the cosmic-ray spectrum 
produced by solar flares is very steep. 

The third problem is to take into account the effect that the refractive 
index n = [1—(@},/w*}!? is not equal to unity for the synchrotron 
mechanisms (the plasma mechanisms are not affected by n ~ 1, if 
W >> Wye). 

The fourth problem is to take into account the exchange of energy 
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between cosmic electrons and cosmic ions. If only heavy ions are injected, 
by low-frequency oscillations, they can in a dense and turbulent plasma 
have the same power-type spectrum if ¢ >> (m,c?) (m,/m,). 

The last problem, which is the most important one for astrophysical 
applications, is to construct a model in which the particle and turbulence 
distributions are spatially dependent. Some of these models can include 
in the centre a very turbulent and dense region in which high-energy 
electrons and ions with power-type spectra are formed. They diffuse into 
the outer regions, excite the turbulent motions, then scatter off them, 
change the distribution function and penetrate to larger distances from 
the central regions. This picture, in principle, can give an answer to the 
spectrum of particles and waves emitted by this region and connect the 
spectrum of particles produced in the inner, optically thick, regions 
with the spectrum of particles leaving the source. One can expect that 
through this diffusion to outer regions the es of the particle distribution 
is lowered and forms a distribution with relatively small ¢,. If this picture 
is constructed, one must face the problem of the energy source of the 
central region, and this, of course, poses other purely astrophysical and 
more general questions. We must emphasise one point. We can say that 
we know that it seems that an acceleration by the induced Compton 
effect can lead to a narrow y-distribution of the particles accelerated 
in the range that corresponds to the values observed. Further development 
of the theory must produce answers to the questions which we have just 
stated. 


Conclusions 


THE preceding text can only represent the present state of the development 
of the theory of plasma turbulence. The author has tried to follow the 
logic, but not the history of this field and, therefore, the references are 
very fragmented and not by any means complete. The essential physical 
statements that the author wants to emphasise finally are: 


1. The plasma properties in the turbulent region are mostly non-linear. 
This raises the possibility of universal plasma properties like a universal 
spectrum that can be independent of the type of instability. 

2. Nevertheless, the turbulence is often weak: W/nT<« 1, and when 
describing the properties of the turbulent oscillation interactions it is 
not possible to expand the non-linear interactions in terms of the turbu- 
lent energy. The elementary excitations such as plasmons and “dressed” 
particles have thus a finite lifetime which is connected with their non- 
linear interactions. 

3. The small low-frequency perturbations in a turbulent plasma have 
quite a different nature because of the frequent turbulent collisions, and 
the dielectric constant that describes such perturbations cannot be 
expanded in terms of the turbulent energy. 

4. The development of a turbulent state is very probable for a plasma 
as a result of the fact that the energy applied has a tendency to disperse 
to the greatest possible degree of freedom. Innumerable numbers of 
different plasma instabilities can bring the plasma to a turbulent state. 
The plasmas in astrophysical conditions must, therefore, often be 
turbulent. This can lead to a way of explaining cosmic-ray origins with 
a universal power-type spectrum. 

5. The development of plasma turbulence can occur as a result of 
development, firstly, of one or a small number of collective modes, with 
a subsequent spread of the energy to other modes by non-linear inter- 
actions as well as by the excitation of many modes at the first stage. For 
the case of the excitation of one mode, the first stage is not turbulent 
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and the turbulence develops as the energy is spread, if the system is 
ergodic. The plasma collective motions seem to be the best test for an 
investigation of the general problems of the development of the random- 
isation process, as well as of the general problems of the possibility of 
a Statistical description of a system. 

Further experimental and theoretical developments of this problem 
seem to be one of the future interests in this field. 
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